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I. ABSTRACT 


An implicit finite-difference scheme is used to solve the laminar 
and turbulent boundary-layer equations for perfect gases and reacting 
gas mixtures in chemical equilibrium. The formulation of the boundary- 
layer equations neglects transverse curvature effects, and the equili- 
brium chemistry model assumes that the element composition across the 
boundary-layer is constant. Thus, injection of a foreign gas at the 
wall boundary cannot be considered. 

The numerical procedure is applied to both internal and external 
flow problems and the results are compared with experimental data and 
other numerical solutions where these data were available. The solu- 
tions for laminar and turbulent flows of perfect gases and laminar 
flow of an equilibrium gas without mass transfer are in good agreement 
with experimental data and/or other numerical solutions. For the case 
of mass transfer at the wall, the numerical solution is in good agree- 
ment with experimental heat-transfer data, but the velocity profiles 
and skin-friction predictions are not in good agreement with the 
available data. Additional experimental data are needed to assess the 
accuracy of the numerical solutions with mass transfer effects. 

The experimental data which were available for turbulent flows of 
an equilibrium gas were in a range of pressure and/or temperature where 
the effects of equilibrium chemistry are not large. With the exception 
of an Integral method of solution for turbulent flow of an equilibrium 
gas in nozzles - which failed to converge for the problem considered - 
other numerical methods were not available for comparison. However, 
the solutions obtained are in good agreement with the limited data 
available. 
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V. INTRODUCTION 


The existing literature on the numerical solution of the laminar 
boundary-layer equations for two-dimensional and axisyraraetric flows is 
extensive. A recent review of the most commonly used techniques for 
solving the laminar boundary-layer equations for non-equilibrium, 
equilibrium, and non-reacting chemistry is given by Blottner [ref. 1], 
Kline et al. [ref. 2] present a similar review of the prediction 
methods used for the solution of the incompressible turbulent boundary- 
layer equations. Examples of the more recent solutions of the tur- 
bulent boundary-layer equations are the reports by Harris [ref. 3] and 
Pletcher [ref. 4]. 

The implicit finite-difference scheme of the Crank-Nlcolson 
[ref. 5] type has been developed extensively by Blottner [refs. 1, 6, 

7] and by Davis [refs. 8, 9, 10] for a wide range of laminar boundary- 
layer flows. This method of solution has been demonstrated to be 
accurate and stable and does not require an excessive amount of com- 
puting time. This type of finite-difference scheme has been used by 
Harris [ref. 3] to solve the turbulent boundary-layer equations for 
non-reacting gases. Harris considered mass transfer at the wall and 
the laminar-turbulent transitional regieme. Cebeci et al. [refs. 11, 
12, 13] used an implicit finite-difference scheme to obtain the solu- 
tion of the turbulent boundary-layer equations. However, the numeri- 
cal procedure used by Cebeci differed considerably from the Crank- 
Nicblson type scheme. The solution of Pletcher used an explicit 
finite-difference calculation procedure based on the DuFort-Frankel 
[ref. 14] scheme. The turbulent solutions of Cebeci and Pletcher also 
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considered only non-reacting chemistry. 

In the references cited above , the authors have considered only 
external flows. Elliott, Bartz, and Silver [ref. 15] have developed 
an integral method of solution for predicting turbulent boundary- 
layer flows in rocket nozzles. Boldman et al. [ref, 16] have applied 
the method to predict turbulent flows in supersonic nozzles. Eden- 
field [ref. 17] has extended the method to predict turbulent flows in 
hypervelocity nozzles in which the gas- was considered to be in chemi- 
cal equilibrium. 

For the lower Mach number perfect ga9 cases considered by Elliott, 
Bartz, and Silver and by Boldman, relatively good agreement between the 
predictions and the experimental data for the heat-transfer distribu- 
tion was obtained. Edenfield found that the method did not predict 
the boundary-layer displacement thickness accurately downstream of the 
nozzle throat and failed to converge for local Mach numbers of about 
16, Thus, for hypersonic nozzle flows, a limit exists for the use of 
the Integral method. Other disadvantages of the Integral method are 
the amount of empirical data needed and the number of adjustable para- 
maters which strongly Influence the results of the predictions. 

As a result of the excellent agreement between experiment and 
theory which has been obtained by Blottner and Davis for both reacting 
and non-reacting chemistry for laminar boundary-layer flows and by 
Harris for non-reacting turbulent flows using the Crank-Nicolson type 
implicit finite-difference scheme, this method of solution was se- 
lected for the present Investigation. Both laminar and turbulent flows 
of perfect gases and mixtures of perfect gases in chemical equilibrium 


2 



are considered for flat plates, wedges, two-dimensional or axi sym- 
metric blunt bodies and nozzles. Hass transfer is considered for the 
case where the injected gas is the same as that of the external flow. 
The primary emphasis has been to obtain solutions for high Hach number 
flows having strongly favorable pressure gradients and highly cooled 


walls 



VI, ANALYSIS 


The equations of notion for laminar or turbulent flow of perfect 
gases and equilibrium gas mixtures are developed in Levy-Lees vari- 
ables and expressed in the general parabolic form necessary for the 
implicit finite-difference solution procedure employed by Blottner 
[ref. 6} and by Davis [ref. 8]. 

Semi-empirical expressions for the turbulent eddy viscosity are 
presented for the cases of a solid wall and a porous wall. The 
finite-difference scheme is developed and the procedures employed for 
determining the initial profile data, and the specification of edge 
conditions are discussed. Definitions of the boundary-layer para- 
maters which have been employed are also given. 

6 . 1 Governing Equations 

The governing equations for laminar or turbulent boundary-layer 
flow of an arbitrary gas in thermodynamic equilibrium or of a perfect 
gas are presented in dimensional variables and transformed to Levy- 
Lees variables. The rate of mass transfer at the wall boundary for 
porous-walls is assumed small in comparison to the boundary-layer mass 
flow and normal gradients are negligible. The boundary-layer thick- 
ness is assumed to be small in comparison to the body radius of curva- 
ture and centrifugal forces are neglected. The coordinate system is 
shown in Figure 1. 

6.1.1 Laminar Boundary-Layer Conservation Equations 

The conservation equations for laminar boundary-layer flows of a 
perfect gas or of a chemically reacting gas mixture in equilibrium are 
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developed In this section* For a multicomponent gas mixture, the lam- 
inar boundary-layer equations are expressed in dimensional variables 
aB; (see Hirschfelder et al, [ref. 18]): 

Continuity 


* * 
3p u r 
* 

3 x 


* * *4 

3p v r 
* 

3 y 


(1) 


where 


j « 0 for two dimensional flow 
j •• 1 for axisymmetric flow 


Momentum 


**3u **3u 

p u + p v £ 

3 x 3 y 



* 

* 3 u 

y 

9 7 


( 2 ) 


Energy 


* * 3 h* , * * 3 h* 

p u + p v * 

3 x 3 y 




3 u 


[3 y 


(3) 


where 


* 

q 



iss 

l 

i«i 


(4) 


and 


* 

j* « £— 
i Pr 


ISS M M 

i -Vi 

J-l 


iii 

iJ 3 y* 


'1 3 1 

* „ * 
r 3 y 


(5) 
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Species 


* * 8 c i * * 3 c i 

p u —i + p v — i 

3 x 3 y 


3 y 


(J*) 


( 6 ) 


Assuming the gas mixture to be in local chemical equilibrium, the spe- 
cies composition is a function of the pressure, temperature, and con- 
centrations of the chemical elements. The diffusion equation for the 
jth chemical element is obtained from equation (6) after multiplica- 
and summing over all species, ISS, as 


tion by 


iM i i 


* * 3 c' 


J 


P u 


+ p V 


* * 3 c^. 


3 x 


3 y 


3 y 


1 J i 

i-1 1 i 1 


(7) 


where 

. ISS . 

c ■ I (8) 

i-1 1 n i x 

Assuming that the element composition, c^ , remains constant across the 
boundary-layer, the conservation of energy, equation (3), can be ex- 
pressed in terms of total enthalpy in the same general form as that 
for a perfect gas. If c^ is constant, the heat transfer can be ex- 
pressed as 


(9) 


where 


K 


* 

* 

* 

3 T 

ja 

1 

i 

s* 

* 



3 y 

* ISS . 

• 

ISS 

M ± M ^ 3 X ± 

a — y h 

y 


rr f i-1 1 

L 

» 

M 2 ^ 3T* 


*1 


X, 

* 


( 10 ) 
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Using the definitions, 


«fc ^ it 

H “ h + “ 


( 11 ) 


h - h (P , T ) 


and the normal momentum equation 


3P 


- 0 


3y 


the heat transfer, equation (9), can be expressed as 


( 12 ) 


(13) 


Pr 


* * 
3H * 3u 
. * - u 
[3y dy J 


(14) 


where 


Pr 


* * 

v 

TT 


C P u 


K 


(15) 


Multiplying equation (2) by u and adding to equation (3) with the use 
of equation (14) gives 


* * 3 H* . * * 3 H* 

p u * + P v £ 

3 x 3 y 


V* 3 H* *f 1 ) * 3 u* 

^ 1 - id u ; 

d V * ft V 


(16) 


3 y LIT 3 y ' a y J 

The conservation of energy as expressed by equation (16) has the 
same form as that for a perfect gas, but for an equilibrium gas the 
thermodynamic and transport properties are determined for the specified 
gas mixture. This approach is limited since the element composition is 
not strictly constant across the boundary layer. This approach cannot 
be used for injection of a foreign gas, since the element composition 
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is not a constant across the boundary layer. Therefore, the case of 
mas 8 transfer Into the boundary layer Is restricted to Injection of 
the same equilibrium gas mixture or perfect gas as that at the outer 
edge of the boundary layer. 

For high Reynolds number flow, the conservation equations are non- 
dimensionallzed by variables, given by Van Dyke [ref. 19], which are of 
order one in the boundary layer. The nondimens ional variables are de- 
fined in the list of symbols. The conservation equations in non-dimen- 
sional variables have the same form as the dimensional equations and 
arei 

Continuity 


3 pur^ , 3 pvr^ 
3x 3y 


0 


(17) 


Momentum 


Energy 


3u . 3u dP e .3 3u] 

pu te + pv 3y " " ~ + a7 3?J 


(18) 


3H . 3H 
pU 3^ + pV 


For a perfect gas. 


3_ [£_ 3H 
3y Pr 3y 


+ V 



(19) 


Pr * constant 
* 

Cp ■ constant 

h - T (20) 

and the viscosity may be expressed by Sutherland's law 
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A more convenient form of the conservation equations for numerical 
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solution is obtained by the introduction of a stream function defined 
as 


Kx.y) - (20 1/2 f(S,n) 


( 25 ) 


where 5 and n are the Levy-Lees transformed coordinates! 


X 

■I 


& 


5<x) " I P e u e >* e r dx 


( 26 ) 


pur 
e e 


n(x,y) - "r - m 
(20 ' 


- I p/p e dy 


( 27 ) 


At the stagnation point of a blunt body the boundary-layer equations 
have a removable singularity. As € * 0, the limiting process gives 


du 2(j+l) 

_e x 




( 28 ) 


n(x,y) - 


. *« +1) d« ? 

, £. £_ dv 

( 25) 1 ' 2 


( 29 ) 


or 


n(0,y) - 


du e 

p e dT <j +1 > 


1/2 


dy 


( 30 ) 


The differential operators expressed in the 5 # n coordinate system 
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are! 


3__ 

3x 


2j 3 , 3 

% "e r 3 a? + \17 


and 


3_ 

3y 


0 U . 'j 8 

( 25) 1/l2 Sn 


Using equations (25) , (31), and (32) gives 


(31) 


(32) 


3 23 

•r J • pv " P u u r J 

6 6 o 


(20 


T 7 ? + ««> l/2 £ ? 


+ (25) 1/2 n x F 


or 


25 f ^ + V + f - 0 


(33) 


where 


V 


2L 


P a u u r 
e e e 


2j 


n f + 
x 




(34) 


At the stagnation point, 



Differentiation of equation (25) with respect to y using equation (32) 
gives ! 


f* ■ u/u 

e 


(36) 
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Evaluation of equation (18) at the outer edge of the boundary- 


layer gives the pressure gradient as: 


dP du 

e m e 

dx P e U e dx 


(37) 


With the use of equations (31) - (37) , the conservation equations 
for laminar boundary layers are; (Blottner [ref. 7]): 

Continuity 

25 F ? + V' + F - 0 (38) 

Momentum 

2?F F + VF* - 8 

Energy 

25F g ? + V g’ - g” +~ g’ + 


1&- F 2 


+ (CF')' 


(39) 



PrJ 


[C’F'F + CF' 2 + CF”] + CFF 1 



(40) 


where 


F - f» 



c - -2iL_ 


p Q u 
e e 


and 


8 


2 £^ 

u e d? 


(41) 


12 



or at the stagnation point 


6.1.3 Turbulent Boundary-Layer Conservation Equations 
In this section, the conservation equations for turbulent flow 
and the semi-empirical formulations of the eddy viscosity model and 
the eddy thermal conductivity are presented. Following the usual 
practice, the symbols H, P, p, u, and v are to be interpreted as time 
averaged properties. The nondimensionalized form of the conservation 
equations are: 

Continuity 


3 pur J , jJJ 
3x 


Momentum 


3u , , T 1 T\ 3u U e ,3 3u 

pu ■?— + (pv + p'v*) — ■ p u —— + -r— u r— - P u'v' 
3x 3y e e dx 3y 3y J 


Energy 


pu u + (pv + p*v r ) 


9H 9 u_ 9H -nrr . f, 1 1 

9? *97 p717-< ,v H + “ ( 1 ■ PrJ 


The solution of equations (42)-(44) requires expressions relating 

the Reynolds shear stress term -u'v’ and v *H ' to the mean variables 

u, v and H, These expressions are obtained by introducing an eddy 
* 

+ * 
viscosity, e , and an eddy thermal conductivity, e^, where it is 

assumed that 


13 


(46) 


and 


or 


where 


-u 


t 


* 

v 


* 

+ 

e 



- e *i£- 

k 3y* 


* + * 

C P e 3H 

■ v H ’ " — * 

t 3y 




(47) 


(48) 


(49) 


The value of the turbulent Prandtl number Is taken to be 0.9 for both 

perfect gases and equilibrium gas mixtures. 

* 

The eddy viscosity* e f is evaluated using the concept of a two- 

* 

•f* 

layer eddy viscosity model consisting of an inner law, e. , valid near 

* 1 

the wall and an outer law, e Q , for the remainder of the boundary layer. 
This procedure has been employed successfully by a number of authors; 
for example, Cebeci, Smith, and Mosinskis [ref. 12], and Harris [ref, 

3]. These authors used expressions for the inner eddy viscosity law 
which were based on Prandtl' s mixing-length concept stated as 



(50) 


where Z is the mixing-length. In the present solution of the turbu- 
lent boundary-layer equations, a number of expressions based on equation 
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(50) have been used, and in addition to these models, an eddy viscos- 
ity based on the Boussinesq [ref, 21] relation 


* 


T 


1 + 




(51) 


has been used for the non-porous wall cases. For the porous wall case 
the eddy viscosity expression considered was based on equation (50) , 
These expressions are given below. 


A, Expressions for the Inner Eddy Viscosity Law 
for no Mass Transfer at the Wall 

The eddy viscosity laws based on equation (50) have been derived 
by analogy with Van Driest *s proposal for the mixing-length. Van 
Driest {ref, 22] considered Stokes' flow for an infinite flat plate 
with periodic oscillations in the plane parallel to the plate to ar- 
rive at an expression for the mixing length. The expression given by 
Van Driest is 

£* * y*[l - exp(y*/A*)J (52) 

or using law of the wall coordinates 

^ y*[l - exp(-y + /A + )3 (53) 


where 


+ 

y 


* 

y 


* 


* 


V 
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* 

u f 


- «>*> 1/2 


By correlation of experimentally determined velocity profiles for 
incompressible turbulent flows in tubes the constants were found to be 

k x - 0.4 

and 


or 



* 

A 


» 26v 


* 

IP J 


- 1/2 


(54) 


Patankar and Spalding [ref. 23] proposed the shear stress at the 

* * 
wall, x , in equation (54) be replaced by the local value x . The 

expression for A using this proposal is 

A* - 26v*(x*/p*)" 1/2 (55) 


The conservation of momentum, equation (44) , for an incompressible 
two-dimensional flow can be expressed as 


3x 

* 

3y 


dP 


ji 

* 


(56) 


dx 


for the region near the wall. Integration of equation (56) and sub- 
stituting into equation (55) gives A* as 
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(57) 


* 

A 


* 


- 26v 



+ 



*' 

* 

dP 


e 

* 

* 

p 

dx 


- 1/2 


or 


A + - 26 [1 - P + yV 1/2 (58) 

where 

* 

+ dP e * , * * 3 , . 

P - 1 v /p (u f r (59) 

dx 


Cebeci and Smith [ref. 11] note that the tern In brackets of equa- 
tions (57) and (58) may become negative for accelerating flows leading 
to the square root of a negative number. To avoid the numerical dif- 
ficulty Cebeci and Smith replaced equations (57) and (58) by 


26\> 


* * 

t * dP 

-2 + 2 £ 

* * * 
p p dx 


1 - 1/2 


(60) 


and 


A + - 26 [ 


|i - pV|3- 1/2 


(61) 


As a result of the difficulties encountered using equations (57), 
(60) and (61), equation (57) was arbitrarily modified by replacing the 
pressure gradient with the absolute value giving 
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(62) 


A ■ 26v 


* 

A 

1 

T * 

dP 


-W + 2_ 

e 


A * 

* 


P P 

> 

dx 



- 1/2 


In a more recent publication, Cebeci [ref. 13] suggested that the 
value of y"*" in equation (58) be replaced by a constant value of 11.8, 
This is an experimentally determined value for the laminar sublayer 
thickness and gives 

A + « 26 [1 - 11.8P + r 1/2 (63) 

Cebeci [ref, 13] does not indicate a procedure for the case of a nega- 
tive square root which occurs for P*" > 1/11.8, 

Relchardt [ref. 24] considered the incompressible continuity equa- 

* 

tion for the fluctuating velocity components to demonstrate that e 
3 1 

varies with y and presented an expression which was obtained by curve 

fitting experimental data of flow in pipes. Reichardt ' s expression for 

the inner eddy viscosity is 




* * 

* *' 

a 

p 0.4 

7 u 

y u f 

^ — xx tan n 

* 


V 

11 V 



k 4 


(64) 


B. Eddy Viscosity Expression for the Case of a Porous Wall 
For a porous wall with pressure gradient, the conservation of 
momentum, equation (44), is approximated for the region near the wall 
as, (Cebeci [ref, 13]): 
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( 65 ) 


* v* * dP* 

dT w * e 

* * T * 

dy v dx 


Solving equation (65) for the shear stress gives 


A 

. dP * 
* e v 


dx v 


T I l“T 7 - 1 + T « e xP T 7 


The damping constant , A , for the porous wall is expressed as 


A + « 26 | - ~ [exp(vpr + ) - 1] + exp(v*y + ) 


Assuming that y at the edge of the laminar sublayer for the case of a 
porous wall Is approximately the same as that for a flat plate without 
mass transfer, Cebeci [ref, 13] used a value of 11,8 for y + and ex- 
pressed A + as 


A + - 26 


■£- [exp (11, 8 v*) 

v 

w 


1] + exp (11. 8 v w ) 


For the case of no mass transfer, equation (68) reduces to equation 
(63), and for a porous flat plate becomes 


A + ■ 26 exp (-5, 9 v*) 


The nondimens lonal fora of the inner eddy viscosity laws for no 


mass transfer at the wall are: 
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Absolute Value of the Pressure Gradient 



( 70 ) 


( 71 ) 


( 72 ) 


( 73 ) 
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Cebeci l~ref. 131 




( 74 ) 


where A is given by equation (63). In equations (70)-(74), X^ and 
are defined as: 


and 


X l" 


2 2 2 2 1 
r J 
e 

TTT 


y p u r 


e VD V (2 ° 


X 


2 


IE 


26p 



(75) 


(76) 


Reichardt [ref, 241 


+ 

e i 




1/2 

* 


1/2' 

c f p 


C f p 

00 

- 4.4 tan h 

y 

*00 


2e VD 

k * 

lly 

t 

2£ VD 

k i 

» 


(77) 


The inner eddy viscosity law for the case of mass injection is (Cebeci 
[ref. 13]) 


+ 

e i 


* 


' 

1 - exp 

_ 

+ 



A 

W > 



|f* 


where A 


is given by equation (68) 


or equation (69). 


(78) 
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Proceeding in the same manner as for the laminar boundary-layer 
equations (Section 6*1.2) , the equations for turbulent flow are ex- 
pressed in Levy-Lees variables as: (Cebeci, Smith* and Moainskis 
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[ref. 12]): 

Continuity 


25 F ? + V* + F - 0 


Momentum 


25 F F 


+ VF 

■n 




+ (C(l + e + ) F') * 


(82) 


(83) 


Energy 


25 Fg 5 + Vg» - f; 


1 + e 


+ Pr_ 
Pr. 


C 

Pr 


1 4* e 


+ Pr. 

Pr. 


u 

g » +-A 

8 H 


1 - [C* FF» 


+ CF’ 2 + CF"] + CFF* [l - 


(84) 


Where 


V - 


2 L 
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) u u r J 
e e e 


£ , + isz* sa r 

(25) 1/2 


( 85 ) 


or at the stagnation point 


V - 


( 86 ) 


pv 4- p’v* 


6,1.5 Boundary Conditions for the Governing Equations 
The boundary conditions at the wall, n ■ 0, and at the outer edge 
of the boundary-layer, n • n e , for equations (34)-(40) and (82)-(86) 
are: 


at n ■ 0: 


0; g 


H h 

-a --a . v - v 

H H * w 

e e 


and 


(87) 


at n ■ n e : F ° 1; g ■ 1 


6,2 Numerical Solution Procedure 

The conservation equations for laminar or turbulent flow have been 
solved using an implicit finite-difference scheme. The numerical method 
is the one employed by Davis [ref. 9] and requires that the governing 
equations be expressed in the general parabolic form 

W" + AjW' + A 2 W + A 3 + A a W ? - 0 (88) 

where W is the dependent variable and the coefficients are functions of 
5, n and W. 

In the following sections, the conservation equations are ex- 
pressed in the general parabolic form, equation (88) , and the finite- 
difference scheme is described. 
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6.2.1 Standard Parabolic Form of the Governing Equations 


The governing equations are expressed in the form of equation 
(88) for a perfect gas or an equilibrium gas mixture asJ 
Momentum 

F" + AjF 1 + A 2 F + A 3 + A 4 F ? - 0 (89) 

where 



A, - - 25F/A 

4 o 


A is defined for laminar and turbulent flow as 
o 

A q ■ C (for laminar flow) 

A ■ C (1 + e + ) (for turbulent flow) 
o 


and 



Energy 


where 


g" + Ajg* + A 2 g + A 3 + A^g ? ■ 0 



(90) 


25 


A 2 - 0 


2 

u , i 

— !=■ i _ i-l 

H r PrJ 


A_ 


f .2 


£L. ff> +Z1 + E1 

CA A A 

o o o 


(for a perfect gas) 


u 

A m C 
A 3 H C 1 
e 


£l + ^i 
c c x 


FF' + F» 2 + FF” 


fk Q (for an equilibrium gas) 


\ - - 2 5F/A 0 


A ■ — (for laminar flow) 
O rr 


A 


o 


C__ 

Pr 


1 + e 


+ Pr 


Pr. 


(for turbulent flow) 



and 


1 


1_ 

Pr 


Continuity 

The solution of the continuity equation is determined by numerical 
integration of the expression 

^e 

V - V - [ (2£ F„ + F) dn (91) 

w j 5 

o 


after each iteration of the momentum and energy equations. The inte- 
gration is performed using the trapizoid rule. 
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6,2.2 Derivation of the Finite-Difference Solution Procedure 
The finite-difference scheme used to solve the boundary-layer 
equations is an implicit method of the Crank-Nicolson type [ref. 5], 
which has been applied successfully by a number of authors; for exam- 
ple, Davis [ref. 9], Blottner [ref. 1] and Harris [ref. 3]. 

The boundary layer is considered as a network of nodal points with 
a varying step size in the normal coordinate direction as shown sche- 
matically in the figure below 




For convenience, equation (88) is expressed at (m,n) and (m+l,n) 


as 


and 



(92) 


n 

m+1 



(93) 


If it is assumed that the dependent variable is known at the 
ppints (m,n) and (m+l,n), Taylor series expansions for gives 


and 


+ <1-X) A? 

m+(l-X) m ' ' (3 5J 


m 


. (1-X) 2 Ag 2 
2 




35 




m 


, 11-1) 3 M 3 
3! 


3 

3 J W 


35 


+• • •+ 


\ * ' 


m 


(94) 


C(l-X) “ Cl " 


f-r 

Na+l 


2 2 
X A5 


3^ 

n 3 3 

X J A5 

f 3 
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7 

3! 

_ _3 

[35 

m+1 



n 


+•• •+ 


m+1 


( 95 ) 


where 



Noting that 



n 

[A] 

n 

f 

3 3 W 

35 2 

0+1 

[h 2 J 

+ AC 
a 

a ? 3 
7 / 


+• • •+ 


m 


and solving equations (94) and (95) for tA , gives 

10+1 


o+l 


VC., « W 11 + A? 

ZQ 


(1 - X) 


3W| n 

3Cj 


+ A 


3W 

nj 


An 


o+l 


. AS 2 
+ 2 


(1 - 2X) 

2 

0 

n 


Lae J 

o 


+ 0(A£ J ) 


(96) 


or using equations (92) and (93), 

Cl “ W o + A5l(1 - X > G - + X 


o 


0+1 J 


+ a(A£ 3 ) (97) 

Davis Iref. 9] approximated aquation (97) as 

Ci - M » + «t» - *> < + >■ Ci 1 ' «» 

and evaluated the A ± in equation (88) at the points (m+l,n). Following 
Davis* formulation, the normal derivatives in equation (88) are re- 
placed by Taylor series expansions for varying step sizes in the n di- 
rection at the point m+l,n as 




a - 2X>i 


'o 
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Ar l 1 ,, 

^ ; m+l An, (An, + An ,) m+1 
n ti n-i 

An - An i An , 

. n n— 1 Tjii _ h Tj n “^ 

An„ An„ , m+1 An, .(An, + An ) m+1 

D n-i n-i n n*ri 

+ 0(An 2 ) (99) 

and 



- - 2 W n ^ — ^ \j a 

An, (An, + An, .) m+1 An, An, , m+1 

n n n-i n n— i 


An ..(An, + An, .) m+1 

n— i n n— i 


W"TJ + 0(An 2 ) 


( 100 ) 


Evaluating the at m+l,n and letting 


and 


W, » W n 
1 in 
n 


W 2 - l W 

n 


gives upon substituting equations (99) and (100) into equation (98) 


A W„ + B W„ + C W„ -D;2<n<N-l 
n 2 , n2 n 2 . , n — — 

n-1 n n+1 


( 101 ) 


where 


2 - A. An 
1 n 

a . n . 

n ” An, , (An, + An, , ) A 
n— i n n— i 
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n 


2 + A (An - An, ,) 
1 n n-l 

n 


AT1 « An r ,-1 

XX XX X 


+ A„ 


X + 


AC 


and 


2 + A, An , 

1 n-1 

Q m S 

n An (An, + An, ,) 
n n n-x 


\ M i 

D - - [W" + A^ W* + A 2 V ] (1 - A) - A^ + -fe-2- 
n n n n n n 


Assuming that 


W 2 " E n W 2 +1 +F n 
n n+1 


(102) 


is valid throughout the boundary-layer (Richtmyer [ref, 26]) ( then 
is given by 

n-1 


W 2 . " Vl W 2 + F n-1 

n-1 n 


(103) 


Using equation (103) in equation (101) and solving for and 

n 

comparing with equation (102) gives 


- C 


n B + A E . 
n n n-1 


(104) 


and 


F - 


D - A F , 
n n n-1 


n B + A E . 
n n n-1 


(105) 
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The values of and are determined by the boundary condition 
for W. 


6.2,3 Spacing of Node Points in the Normal Coordinate Direction 
The finite-difference solution procedure has been developed for a 
variable spacing of the node points in the normal (n) coordinate direc- 
tion. This permits a close spacing of points in the region near the 
wall where the variation of fluid and dynamic properties is greatest. 
The procedure employed is that given by Cebeci, Smith, and Mosinskls 
[ref. 12], 

Using this procedure, the ratio of the adjacent Intervals is a 
constant expressed as 


k 


An 


An 


n-1 


(106) 


The distance to the nth point measured from the wall boundary is 
given by 


An. 



(107) 


where 


An, 



(108) 


N is the number of strips in the boundary-layer, and n g is the location 
of the boundary-layer outer edge. 

The values of the constant k for laminar or turbulent flow were 
determined by numerical experiments, A number of solutions were ob- 
tained for laminar and turbulent boundary layer flows at supersonic 
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conditions considering both adiabatic and cold wall cases. The value 

of k was varied over the range from 1 to 1.5 with n e varied from 4 to 

12 for laminar flow and from 75 to 200 for turbulent flow. The results 

of these test cases were compared with experimental data to determine 

a value of k which gave good agreement with experimental data and 

where the solution showed little change for different values of r> . 

e 

The values selected were: 

k “ 1.04 and n ■ 6 for laminar flow 
e 

k m 1,09 and n - 100 for turbulent flow 

.... e 

The above values correspond to a value of N - 100 in equation (108). 

Similar tests were made varying N from 50 to 500, An N of 50 was 
found to be unsatisfactory for most cases, but values of N greater than 
100 did not improve the solutions which were obtained, 

6,2.4 Convergence Criteria 

A suitable convergence test may be established for laminar 
boundary-layer flows by comparing F^ at successive iterates of the 
solution. This type of convergence test was employed by Davis [ref. 

9] and Cebeci and Smith [ref, 11]. 

For turbulent flows, Cebeci and Smith found that a second require- 
ment based on successive iterates of the boundary-layer displacement 
thickness was necessary to obtain satisfactory solutions. They con- 
tinued the iteration procedure until both tests were satisfied. 

For the present calculations provisions were made to establish 
convergence by comparison of successive iterates of both F^ and g^ or 
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by comparing successive iterates of both F and g at all points in the 
boundary layer* The convergence test are expressed as 
Derivative Test 


and 


1 - 


„ , (k+1) 



F ,(k) 

F W 


± e d 


1 - 


, » (k+1) 

^ ; 

..ao 

% 


± E d 


Function Test 


and 


(k+1) 


1 - 


(k) 


± e f 


n 


1 - 


(k+1) 

!7k r 


i e f 


where k is the iteration number and and are prescribed values. 
The function test with * 0.01 and the additional requirement that 
the number of iterations at the new location be two or greater was 
found to be satisfactory for both laminar and turbulent flows, A 
value of « 0,001 was used with the derivative test, but this test 
is in general not recommended for turbulent boundary-layer calcula- 
tions. 
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6.3 Specification of Body Geometry 


The geometry of a given configuration is considered as a series 
of segments of fourth order and requires that the coordinates, z, x, 
and r be given in tabular form, Since transverse-curvature effects 
have been neglected in the governing equations, provisions have not 
been made to distinguish between different body sections. Thus, some 
inaccuracies are introduced in regions near the intersection points of 
different contours. For problems requiring a specified pressure 
and/or a wall temperature distribution, these data are entered in tab- 
ular form at the same points as the geometry data. The temperature is 
entered in °R and the pressure is given as P/P Q or P/P^. 

To insure accurate interpolation and/or differentiation of the 
temperature and pressure data in regions of large gradients, it is nec- 
essary to have a close spacing of the coordinate data. 

For nozzles, at least 30 points should be entered in the throat 
region, and for blunt bodies, at least 50 points should be entered 
for the nose section. A maximum of 500 points may be tabulated, and 
the minimum number permitted is 5. Best results were obtained when 
the maximum ratio of adjacent step sizes in z were not greater than 
1.25. 

The body shapes which have been considered are blunt bodies, 
wedges, flat plates, and nozzles. 

6.4 Fluid Properties at the Outer Edge of the Boundary-Layer 

In the present method for solving the boundary-layer equations, 
the effects of mass entrainment of the outer inviscid vortical stream- 
line into the expanding boundary-layer on a blunt body is neglected. 
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Therefore, the edge conditions may be specified by conditions on the 
body surface for the inviscld flow of the gas. The procedures em- 
ployed for specifying the outer edge conditions for the different geom- 
etries and gas models is discussed in the following sections, 

6,4,1 Axisymmetric Nozzles 

The edge conditions for nozzle flows of a perfect or an equlli-* 
brium gas may be specified by assuming a one-dimensional expansion or 
the expansion can be determined for a specified pressure distribution 
where the pressure is given in tabular form as indicated in Section 
6.3, The procedures employed for the different cases are discussed 
below, 

A, One-Dimensional Expansion of a Perfect Gas 
For a one-dimensional expansion of a perfect gas, tables of Mach 
numbers and the corresponding area ratios are computed within the 
boundary-layer computet program using the relation 
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at intervals of 0,01 .in Mach number for M < 2 and in intervals of 0,05 
for M ^ 2, 

After the above tables have been generated, the Mach number at 
the nozzle exit is determined by five point interpolation in the tables 
of area ratio and Mach number with the area ratio as the independent 
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variable. The exit Mach number is denoted by M^, 


The free-stream pressure, temperature, velocity, density, vis- 
cosity, and Reynolds number are then computed using the relations 
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where P re £ » ' r rQ £ and L are reference values of viscosity, tempera- 
ture and length. 

Corresponding to the points in the geometry arrays, tables of the 
edge pressure and velocity are computed in nondimensionalized form 
using the expressions 
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and is the local Mach number determined for the area ratio at the 
given point in the geometry arrays. 

At a local solution station x ** x^ + Ax, u g and are deter- 
mined by interpolation, and T^, and are computed using 
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du 

The velocity derivative is evaluated numerically. 
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B. Pressure Distribution Specified for a Perfect Gas Solution 
For solutions where the pressure distribution is input to the 
computer program, the outer edge conditions are computed in the same 
manner as for the one-dimensional expansion solution with the exception 
that the Mach number table is computed from the input pressure distri- 
bution using 
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After computing the Mach number array using equation (125) , equa- 
tions (116)- (124) are evaluated using the local Mach number data cor- 

du 

6 

responding to the given pressure distribution, and is evaluated 
numerically. 


C. One-Dimensional Expansion of a Reacting 
Gas in Chemical Equilibrium 

' The boundary-layer edge conditions for the equilibrium gas case 
were determined by the solution of the inviscid equations of motion 
for the given mixture of perfect gases. The data which must be given 
for the boundary-layer solution are: 

A/A , M, 0*, ?, and ft 

U 

These data are tabulated for use with a table-look-up procedure. For 
a one-dimensional expansion of the gaa mixture, the area ratio is used 
as the independent variable to determine the local edge conditions for 
the given body geometry. Since the area ratio is not a monotone func- 
tion over the length of the nozzle, it was found to be more satis- 
factory to create a secondary expansion data set on a scratch unit ex- 
pressed in terms of z or x. The locations of z and x were determined 
for a given A/A £ by interpolation and the secondary data set was writ- 
ten as unformatted records in the order 

* '\j 

z, x, M, U , P and h 

This operation is performed within the boundary-layer program and re- 
quires no additional preparation of the expansion data. In addition 
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to the expansion data tape, a tape unit Is required for the tables of 
thermodynamic and transport properties of the equilibrium gas mixture. 
These tables are written at constant $ with the temperature decreasing, 
The data in each table are written as unformatted records in the form 
(Pj T*, ft, p, p*, 0 ^, Pr). The expansion data and the tables of 
thermodynamic and transport properties were obtained using two modified 
versions of the computer program developed by Lordi et al. [ref, 20], 

A description of these modifications are given in Appendix A, 

After the secondary expansion data set has been written, local 

it. f\, f\j 

values of M, U , P and h are determined by interpolation with either 
z or x as the independent variable. With the local values of ^ and h 

* *\j * 

as independent variables, the local values of T , p, p , c^ and Pr are 

found by interpolation in the tables of thermodynamic and transport 

properties. The nozzle exit conditions are taken as the free-stream 

conditions and the reference conditions retain the same definitions in 

* 

the dimensional form with the exception of b re £» The reference vis- 
cosity for the equilibrium cases may be chosen arbitrarily and does 
not necessarily correspond to the reference temperature and reference 
pressure. The reference viscosity employed is computed using Suther- 
land’s law. 

D, Equilibrium Gas Solutions with the Pressure Distribution Given 

A/ 

If the pressure distribution is given, the local value of P is 
determined by interpolation in the given pressure distribution table. 
With this value of $ as the Independent variable, local values of M, 

U and h are determined by interpolation of the data on the expansion 
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data tape. It is not necessary for P to be a monotone function. The 
remaining edge conditions are determined by interpolation in the 
tables of thermodynamic and transport properties with P and h as the 
independent variables as in paragraph C above. 

6.4.2 Blunt Bodies 

The solution of the boundary-layer equations for flows over blunt 
bodies requires that the pressure distribution be specified. The re- 
maining edge conditions are determined for a perfect gas or an equili- 
brium mixture of perfect gases as discussed below. 

A. Isentropic Expansion of a Perfect Gas Along 
the Body Streamline From the Stagnation Point 
For these solutions, the pressure distribution is entered in the 
table as P/P^ where P^ is the stagnation pressure behind a normal 
shock expressed in nondimensionalized form as 


P' 
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(y+1) M l 

Y-l 
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2y M I - (y-1) 
00 ■— 1 



The reference conditions are based on the flow properties ahead of the 
bow shock. Since the inviscid flow along the body streamline is isen- 
tropic, the edge conditions are computed in the same way as for a noz- 
zle with the pressure distribution given and is described in paragraph 
C of Section 6.4,1 above. It is noted that the expansion is from the 
stagnation conditions behind the normal shock. 
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B» Isentropic Expansion of an Equilibrium Gas 
Alone the Body Streamline from the Stagnation Point 
For this case, the equilibrium expansion data are determined for 
the stagnation pressure and temperature behind a normal shock for the 
given mixture of perfect gases and free-stream conditions. The free- 
stream and reference conditions with the exception of the reference 
viscosity are determined by the fluid properties ahead of the bow 
shock. The edge conditions at local points along the body are deter- 
mined in the same manner as for an equilibrium solution for a nozzle 
with the pressure distribution given and is discussed in 6,4,1 D, 

6,4,3 Wedges and Flat Plates 

The edge conditions for wedges and flat plates are constants cor- 
responding to conditions behind an oblique shock or the free-stream 
conditions respectively. The computer program for the solution of the 
boundary-layer equations is suitable for both wedges and flat plates 
if a perfect gas is considered. The solution of wedges for an equili- 
brium gas requires modification of the boundary-layer computer program. 
Thus, only flat plates are considered for equilibrium chemistry solu- 
tions, The procedure for determining the edge conditions for a per- 
fect gas and an equilibrium gas mixture is discussed below. 

A, Edge Conditions for Flat Plates and 
Wedges for a Perfect Gas Solution 
The pressure, temperature, and velocity at the outer edge of the 
boundary-layer are computed from the oblique shock relations and are 
expressed in nondimens ional form by the relations 
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The density and viscosity are computed using equations (122) and (123) 
or (124). 


B. Flat-Plate Solutions for Equilibrium Gases 
The solution of flat plates in an equilibrium gas requires the 
use of the tables of thermodynamic and transport properties for the 
given gas mixture (Appendix A) . The free-stream pressure, temperature, 
and velocity are input data to the boundary-layer computer program. 
Using the free-stream pressure and temperature as independent vari- 

•ft ft ft ft 

ables, the values of h , p . u , and c are determined by the use of 

* CO* ^CO* p J 

"CO 

interpolation in the tables of thermodynamic and transport properties. 


6,4,4 Evaluation of the Longitudinal Coordinate-^ 

The coordinate £ is evaluated by equation (26) numerically using 
a modified Simpson's rule (Davis [ref. 9]), and is expressed as 


44 



(129) 


^m+1 


+ 

m 


tu 


m 




+ 4(U 

e m+l/2 


2j 

r m+l/2 P < 


) 


! mfl/2 m+1/2 


+ U 

e m+l 



p v ] Ax/6 
m+1 urfl 


where the subscripts refer to the edge conditions at the points x^, 

Ax 

x^ + Y“ t and x m + Ax respectively, 

6, '5 Solution for the Initial Profile Data 
The equations governing the laminar or turbulent boundary-layer 
flow of a perfect gas or an equilibrium mixture of perfect gases 
reduce to a system of ordinary non-linear differential equations at 
x * 0* However, for fully developed turbulent flow, the eddy viscosity 
cannot be evaluated at x ■ 0, for flat plates or nozzles. The limit- 
ing form of the differential equations are employed at x ** 0 for blunt 
body flows. The procedures used to obtain the starting profile data 
are discussed below. 


6,5,1 Initial Profiles for the Solution 
of the Laminar Boundary-Layer Equations 
To start the solution of the boundary-layer equations, initial 
guesses of the profiles for the dependent variables F, g, and V are 
required. At the leading edge or stagnation point, the equations are 
ordinary differential equations, and the initial profile data are 
determined by an iteration procedure using the implicit finite- 
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difference scheme developed in 6.2* 2 . 

Initial guesses are made for the profiles and for the first itera- 
tion* These estimates are denoted by the subscripts ( and ( > c 
respectively and the k+1 iterate is denoted by the subscript ( ) j • 

The initial guesses of the profile data are assumed to have the 
forms given below j 


at n < n 


at n ■ n 
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The solution of the continuity equation is assumed as 


V - 

c 
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and the temperature distribution is computed from the g profile using 
the definitions 
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Solving for 0 gives 
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For a perfect gas 
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and the Chapman-Rubesin factor, C, is assumed unity across the boundary- 


layer. 

Using the above guesses of the profile data, the governing equa- 
tions are solved to determine F 2 , F^, FJJ, g 2 » g 2 , and g' 2 . The sub- 
scripted variables ( and ( ) c are set equal to the variables with 
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r 0 

subscript ( )^ and new values of 0, 0', — , V c , and C are then 

computed where 
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C - -££— = e 1 ^ 2 - l 1 .- * ., ? ! (Sutherland’s law) 
p e P e (0 + c) 

or 

C = g (w-1) (power law) 


The procedure is repeated until the convergence criteria of section 
6.2.4 is satisfied. 

For an equilibrium gas, the initial guesses of the profile data 
are assumed to be the profile data corresponding to the converged 
solution for a perfect gas described above. Using the perfect gas 
profiles of F and g, the remaining thermodynamic and transport prop- 
erties are determined by use of the table-look-up procedure. Deri- 
vatives of 0 and C are evaluated numerically for the equilibrium gas. 
The iteration procedure is the same as for the perfect gas case, 

6,5,2 Initial Profiles for the 
Turbulent Boundary-Layer Equations 
The solution of flat plates, wedges, or nozzle flows assuming 
fully developed turbulent flow assumes that the profiles at x * 0 and 
x “ 0.001 are similar. This procedure has been adopted because the 
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value9 of e + are zero at x « 0, With this exception and assuming an 

"I* 

initial profile for e as zero, the starting procedure is the same as 
for the laminar boundary layer. 

An alternate method for starting these solutions is to use the 
laminar starting profile at x ■ 0 and assume an instantaneous transi- 
tion to turbulent flow at x > 0,001, Both methods have been used and 
the differences in the resulting solutions are insignificant. 

For blunt body solutions, the turbulent starting profiles are 
determined at x » 0 using the limiting form of the governing equations. 

The initial guesses of the starting profiles for the equilibrium 
solution is determined from the converged profiles for a perfect gas. 
Using these profiles, the iteration procedure is continued using the 
equilibrium gas properties until the convergence tests are satisfied, 

6,6 Boundary-Layer Paramaters 

The definitions of the boundary-layer thicknesses, skin-friction 
coefficient, heat transfer, heat transfer coefficients, and Stanton 
numbers which have been used in the boundary-layer calculations are 
presented in this section. Both the dimensional and nondimenslonal 
expressions are given. 

Velocity or Boundary-Layer Thickness 

* 

The velocity thickness, 6, is assumed to be the value of y at 

* 

■ 0.995, and is determined by interpolation in the velocity profile 
u 

array. 


Incompressible Displacement Thickness 
The incompressible displacement thickness is computed using the 
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two-dimensional definition 




Compressible Displacement Thickness 


Two-dimensional 
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Axisymmetric : 

The axisymmetric compressible displacement thickness is approxi- 
mated as; (Cebeci and Mos inskis "[ref. 27]): 
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Momentum Thickness: 
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Heat Transfer Rate: 

In dimensional variables 
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In nondimens ional form 
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Stanton Number Definitions: 

A. Based on free-stream conditions 
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B. Based on edge conditions 
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Heat Transfer Coefficients: 

A. Based on free-stream conditions 
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Skin- Friction Coefficients: 

A. Based on free-stream conditions 
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B. Based on edge conditions: 
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VII. RESULTS AND DISCUSSION 


The implicit finite-difference scheme developed in Chapter VI has 
been applied to obtain solutions of laminar and turbulent boundary- 
layer flows of perfect gases and equilibrium gas mixtures over flat 
plates, an hyperboloid, a spherically blunted cone, and In axisym- 
metric nozzles. 

Solutions of perfect gas turbulent flows using different expres- 
sions for the inner eddy viscosity law are presented in Section 7.1, 
and are compared with experimental data and/or other numerical solu- 
tions for cases where these data were available. Solutions using the 
equilibrium gas option are presented in Section 7.2, and the solution 
of boundary-layer flows with mass injection are presented In Section 
7.3. 

The numerical solutions presented assume either fully developed 
turbulent or -laminar flow; however the boundary-layer computer pro- 
gram described by Miner, Anderson, and Lewis [ref, 28] provides 
options for either an instantaneous or a continuous transition from 
laminar to turbulent flow. The continuous transition model is based 
on the experimental results of Owen [ref. 29] and is discussed in the 
report on the computer program. 

7.1 Perfect Gas Turbulent Boundary-Layer Flows 

Solutions of turbulent flows over flat plates and in axisymmetrlc 
nozzles using the inner eddy viscosity laws of Van Driest, equation 
(70), Cebeci and Smith, equation (72), absolute value of the pressure 
gradient, equation (73), and Reichardt, equation (77), are presented 
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in this section. The inner eddy viscosity law proposed by Patankar 
and Spalding, equation (71), was found to be unsatisfactory for flows 
having a significant favorable pressure gradient and is not considered. 

7,1,1 Perfect Gas Solutions for Turbulent Flows Over Flat Plates 
For flat-plate flows the inner eddy viscosity laws expressed by 
equations (70)-(74) are Identical and are referred to as the Van 
Driest model. The inner eddy viscosity law given by equation (77) is 
referred to as the Relchardt model. 

Three flat-plate solutions are presented corresponding to case 
numbers 20, 26, and 62 of the experimental data given by Coles [ref. 
30], The results of the present numerical method of solution are 
compared with Coles' experimental data and the solutions of Van Driest 
[ref. 22] and Dorrance [ref, 31], The free-stream conditions for the 
three cases are given in Table I. The free-stream Mach numbers for 
cases 20, 26, and 62 were 3,701, 2,578, and 4.544 respectively, and 
the plates were assumed to be adiabatic. 

Boundary-layer displacement thicknesses predicted using the Van 
Driest and Relchardt eddy viscosity models are compared with Coles' 
experimental data in Figure 2, The displacement thickness predicted 
using the Van Driest and Relchardt eddy viscosity model were found to 
differ by less than 3% for the three cases, but the predicted values 
of displacement thickness differ from Coles' experimental data by as 
much as 20%, The numerically determined displacement thickness is less 
than the experimental value in all three cases. The numerical result 
is in good agreement with experiment for Case 62. 
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The numerically determined velocity profiles at x ■ 21*48 in. for 
Case 20 are presented in Figures 3 and 4, The results using either 
eddy viscosity model are in good agreement with Coles’ experimental 
data. Figure 4 shows the solution given by Van Driest [ref, 22], 

Figure 5 shows the eddy viscosity profiles predicted using the Van 
Driest and Reichardt eddy viscosity models at x * 21,48, The dif- 
ference in the two models is not significant. 

The predicted skin-friction coefficients for the three cases are 
compared with the solution of Dorrance [ref, 31] and Coles’ experi- 
mental data in Figure 6, The skin-friction predictions for Cases 20 
and 26 are in excellent agreement with both the experimental data and 
Dorrance 's solutions. For Case 62, the present method of solution is 
in good agreement with the experimental data points for Reynolds num- 
bers greater than 3 x 10^, but for Re ro - 1,7 x 10^ the present numer- 
ical solution predicts a skin-friction coefficient which is approxi- 
mately 15% lower than the experimental value. The present method of 
solution is in better agreement with the experimental data for this 
case than the solution of Dorrance, 

The above results are representative examples of turbulent flat- 
plate solutions and in all cases considered the numerical results were 
not significantly influenced by the choice of the inner eddy viscosity 
law. However, the use of the Reichardt inner law reduced the comput- 
ing time by a factor of approximately 10 (see Section 7 , 1 , 3 ), 

7,1,2 Turbulent Flow of Perfect Gases in Axisvmmetric Nozzles 
The present method of solution has been used to solve three axisym- 
metric nozzles, and the results are compared with the integral method 
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of solution developed by Elliott, Bartz, and Silver [ref, 15] and ex- 
perimental data where the data were available. The inner eddy vis- 
cosity laws of Van Driest, equation (70), Cebeci and Smith, equation 
(72), absolute value of the pressure gradient, equation (73), and 
Reichardt, equation (77), have been considered for the sample case 
given by Elliott, Bartz, and Silver [ref. 15], 


A, Elliott, Bartz. and Silver Sample Case 
The problem considered consists of a 30° conical inlet section, a 
circular arc throat section with a throat radius of 0,885 in, and a 
15° conical divergent section. The reservoir pressure and temperature 
were 300 psia and 4500°R respectively. The specific heat ratio, y» 
was assumed to be 1,2 and the wall temperature was assumed to be a 
constant value of 1145°R, The data for this case are given in Table 
II, and the nozzle geometry is shown in Figure 7, 

Eddy viscosity profiles predicted at the nozzle throat using the 
four eddy viscosity models are shown in Figure 8, In the region near 
the wall, n < 0,4, the eddy viscosity predicted using the four inner 
laws are essentially identical. The Cebeci and Smith inner law (72) 
predicts a second zero value of e + at n - 1,2 resulting from the pres- 
sure gradient term y/p being equal and opposite in sign to T w /p . 

It is noted that the solution using the Patankar and Spalding expres- 
sion, equation (71), fails upstream of the throat as a result of the 
square root term being negative. By replacing the sum 



P 


dP 
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by 






p dx 


Cebeci and Smith [ref. 11] were successful in obtaining solutions to 
problems having a significant pressure gradient. This model* equation 

(72) * results in substantially smaller values of e + for most of the' 
boundary-layer than is given by the Van Driest Law, equation (70), the 
absolute value of the pressure gradient equation (73) and the Reichardt 
Law, equation (77). The values of e + corresponding to the Van Driest 
expression and the expression using the absolute value of the pressure 
gradient differ by less than 3%, and the e corresponding to the outer 
law, equation (79), n > 1.7, also differs by less than 3%, In the 
sublayer, the Reichardt expression, equation (77), predicts values of 

e which are 10 to 20% larger than those given by equations (70) and 

(73) , but above the match point with the outer law, n » 1,4, e is 
nearly the same as when using equation (70) or equation (73) for the 
inner law. It is emphasized that the outer eddy viscosity law is the 
same expression in all cases. 

The resulting profiles of velocity and temperature differed by 

less than 5% for the solutions obtained using the different viscosity 

* 

laws. The boundary-layer paramaters 6,6,0, c^, etc. also agreed 

CO 

to within 5%. If the inner law of Cebecl and Smith is excluded in the 
comparison, the resulting solutions differed by less than 3%. These 
results indicate that including the pressure gradient term in the inner 
law had little influence upon the eddy viscosity profiles, and had 
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essentially no influence upon the resulting solution. 

As a result of the above, and the formulations of equations (71) 
and (72) for the inner law, these expressions were dropped from further 
consideration. The results which are discussed below were obtained 
using the Van Driest or Reichardt expression for the inner eddy vis- 
cosity law. However, as in the previous cases, use of the Reichardt 
inner law resulted in a substantial reduction of the computer time re- 
quirement (see Section 7,1.3). Since the differences in these solu- 
tions were insignificant, these data are shown as a single curve. 

The heat transfer coefficient, boundary-layer, momentum, and dis- 
placement thicknesses predicted by the present numerical solutions are 
compared in Figures 9-12 with the results obtained by Elliott, Bartz, 
and Silver [ref, 15] using an integral method of solution. The solu- 
tions using the present numerical method and the integral method differ 
by up to 30%; however it should be noted, that the solutions using the 
integral method can be varied over a wide range by changing the assump- 
tions which must be made relating to the nominal entrance conditions. 
Since the starting profiles for the present method are determined from 
the solution of the governing equations, direct comparison of the two 
methods is not possible, also it is not clear what assumptions should 
be made for the initial conditions necessary for the integral method 
if experimental data were not available before the calculations are 
made. Because of the arbitrariness of the solution using the Integral 
method, the agreement of the two solution procedures is considered, to 
be adequate. 
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B. NASA-Lewis Case 2a (Boldman et al. fref. 16 P 


The nozzle configuration consist of a 30° conical convergent sec- 
tion, a circular arc throat section, and a 15° conical divergent sec- 
tion. The throat radius is 0.746 in., and the stagnation pressure and 
temperature were 300 psia and 970°R respectively. The test conditions 
are given in Table III. The geometry is shown in Figure 13, and the 
pressure distribution corresponding to a one-dimensional expansion and 
the experimental data are shown in Figure 14. The experimentally 
determined wall temperature distribution is shown in Figure 15. 

Solutions to this problem were obtained using the Van Driest and 
the Reichardt expressions, equations (70) and (77), for the inner eddy 
viscosity laws, and as in the previously discussed Elliott, Bartz, and 
Silver sample case, the differences in the solutions were insignifi- 
cant. However, the solution using the Van Driest inner law required 
approximately 10 times more computing time than was necessary using 
the Reichardt inner law (see Section 7.1.3), 

The predicted heat-transfer coefficient using the present method 
of solution is compared with the experimental data and the solutions 
obtained using the Elliott, Bartz, and Silver integral method in Figure 
16, The present method of solution is in excellent agreement with the 
experimental data in the throat region and downstream. Differences of 
up to 20% between the predicted and experimentally determined heat- 
transfer coefficient are noted in the subsonic region of the nozzle. 

The near discontinuous change in the experimental value of the heat- 
transfer coefficient at z = 1.97 is the result of the temperature tabu- 
lated in this region (see Figure 15), and also the experimental pres- 
sure data were not smooth in this region. For the present calculations. 
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these data were smoothed in the region 1.9 <_ z <, 2.5. The integral 
method of solution is seen to reflect a strong dependence upon the 
starting condition assumptions. The two solutions presented using the 
integral method differ from each other by as much as 50%. 

C. AEDC Hotshot Wind Tunnel Nozzle 

The problem considered corresponds to the case referred to as 
"Hotshot 1" by Edenfield [ref, 17]. The nozzle geometry, wall enthalpy 
and pressure distribution are shown in Figures 17-19. The operating 
conditions are given in Table IV. In this section the experimental 
pressure distribution shown in Figure 19 was not used. The results 
presented assume a one-dimensional expansion of nitrogen and the speci- 
fic heat ratio, y, was taken as 1.4. The experimental wall enthalpy 
distribution, Figure 18, for 20 milliseconds of tunnel operation was 
used. The stagnation pressure and temperature were 11,500 psi and 
5400 °R. The nozzle consist of 10° and 5 e converging conical inlet 
sections, a faired circular arc throat section, and a 5° conical di- 
vergent section. The throat radius is 0.055 in. 

Attempts to solve this problem using the Van Driest expression, 

equation (70) , for the inner eddy viscosity resulted in an instability 

at z/r ~ 100.* When the Reichardt inner eddy viscosity law, equation 

(77), was used a complete solution was obtained. The solutions ob- 
it 

tained are compared for values of z/r up to 100 in Figures 20-23. 

£ 

Figure 20 shows the eddy viscosity profiles for z/r in the range 
from 7 to 80, Differences In the eddy viscosity using the two models 

*This applies only if the velocity derivative in equation (70) is 
evaluated implicitly (see Section 7.1.3). 
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differ from 5 Co 20% with the larger differences being aC stations near 
the nozzle throat. Heat-transfer rates predicted using the two eddy 
viscosity models are shown in Figure 21, For this problem, the Van 
Driest inner law resulted in heat-transfer rates which are approximately 
10% higher than the predicted heating rates using the Reichardt inner 
law in the region near and upstream of the throat. Downstream of the 
throat, the Van Driest model oscillated about the solution obtained 
using the Reichardt inner eddy viscosity law. The displacement thick- 
ness predictions using the two inner eddy viscosity laws were nearly 
identical and are shown in Figure 22, 

The computer time required to obtain a solution to a given axial 
. * 

position, z/r , is shown in Figure 23. The curves illustrate the time 

requirement in terms of the total number of iterations. The use of 

the Van Driest inner law required approximately five times as many 

■k 

iterations to obtain a solution to z/r ■ 100 as was necessary for the 
complete solution using the Reichardt inner law,* 

The problems discussed in this section have demonstrated that the 
choice of the inner eddy viscosity law had little influence upon the 
resulting solution for cases where convergence was obtained, and that 
the computing time requirement was considerably less when the Reichardt 
expression was used for the inner law (see Section 7,1,3), For the 
cases where experimental data were available, the results obtained with 
the present method of solution were in good agreement with these data, 

*This applies only if the velocity derivative in equation (70) is 
evaluated implicitly (see Section 7,1,3). 
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7.1.3 Revised Evaluation of the Inner Eddy Viscosity Laws 
After the solutions discussed in Section 7.1 and 7.2 were ob- 
tained, the more recent expression for the inner eddy viscosity law 
given by Cebeci [ref. 13] became available. This relation, equation 
(74) , was used in the boundary-layer computer program with and without 
including the pressure gradient term; equation (74) is the same as 
that given by Van Driest, equation (70), when the pressure gradient 
term is neglected. The results obtained using equation (74) with the 

pressure gradient term included were essentially the same as that given 

* 

by the Van Driest expression and resulted in instability at z/r - 100 
for the AEDC nozzle. 

In all of the results which have been presented up to this point 
using Van Driest's and Cebeci' s expressions for the inner eddy vis- 
cosity, the velocity derivative was evaluated at the station x + Ax. 

It was noted that in the cases which resulted in converged solutions, 
the solutions obtained were relatively unaffected by the inner eddy 


viscosity law which was used, and that a large number of iterations 
were necessary if viscosity laws based on Van Driest's [ref. 22] modi- 


fication of Prandtl's mixing length were used. The difficulty appeared 


to be in the evaluation of 


>F~ 1 


To determine the influence of the 


Implicit evaluation of the velocity derivative, the approximation 



was made. With this replacement, the Van Driest inner law, equation 
(70), resulted in solutions which were essentially identical to the 
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solutions obtained using the Reichardt expression, equation (77), and 
the computer time requirement was nearly the same for both eddy vis- 
cosity laws. When the pressure gradient was included in equation (74), 
the results obtained showed that the pressure gradient term had little 
influence upon the solution. For the cases considered, the use of 
equations (70), (74), and (77) resulted in solutions which agreed to 
within + 3%. Therefore, the results comparing the eddy viscosity laws 

presented in Section 7,1 should be considered as the influence of the 

3F 

implicit evaluation of — rather than the differences in the eddy vis- 
cosity expressions. 

The results presented in' Section 7.2 have been obtained using the 
Reichardt expression for the inner eddy viscosity, and as mentioned 
above, the solutions are not significantly affected by the choice of 
the inner eddy viscosity expression. Equations (71) and (72) were not 
considered in this evaluation. 

7,2 Laminar and Turbulent Boundary-Layer Flows of Reacting Gas 
Mixtures in Chemical Equilibrium 

For the perfect gas solutions presented in Section 7.1, it was 
possible to obtain experimental data for most of the problems con- 
sidered. However, it has not been possible to obtain experimental data 
in the range of pressure and/or temperature where the departure from 
the perfect gas condition is significant. For the problems considered 
in this section, both the perfect and the equilibrium gas solutions 
are given for the cases where significant differences were found to 
exist. 

To check the accuracy of the table-look-up procedure, a number of 
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problems were solved at conditions corresponding to the range of pres- 
sure and temperature where the perfect gas relations are correct. For 
the cases considered the solutions with the table-look-up and the per- 
fect gas options agreed to within +3 %. 

7.2.1 Equilibrium Gas Solutions for 
Laminar and Turbulent Flows Over Flat Plates 

Two flat plate cases were considered corresponding to the experi- 
mental data given by Hironimus [ref. 32]. The free-stream Mach number 
and temperature for Case 1 were 7.391 and 533°R and for Case 2, were 
7.58 and 339°R respectively. The wall temperature was 533°R for Case 
1 and 355°R for Case 2. The test conditions are given in Table V. 

Plots of Stanton number based on edge conditions are shown for both 
fully developed turbulent and laminar flow for Case 1 and Case 2 in 
Figures 24 and 25. The solutions shown were obtained with the table- 
look-up procedure. The perfect gas solutions for these cases differed 
from the equilibrium solution by less than 3% and are not presented. 

Downstream of the transition from laminar to turbulent flow the 
agreement of the predicted and the experimentally determined Stanton 
number data is within 10% , and in the region of laminar flow the agree- 
ment between the numerical and experimental data is good. 

7.2,2 Equilibrium and Perfect Gas Solutions for 
Laminar Flow Over a Hyperboloid 

A 10° half-angle hyperboloid at an altitude of 100,000 ft. with 
a wall temperature of 1400°K was considered for a free-stream Mach 
number and temperature of 20.178 and 226,98°K, This case corresponds 
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to Case A of the AGARD test cases (Lewis [ref. 33]). The flight con- 
ditions are given in Table VI. The body geometry and pressure distri- 
bution are shown in Figures 26 and 27. 

The solutions of both perfect air and equilibrium air have been 
obtained for the laminar esse only. The perfect gas solution is com- 
pared with the results of Adams [ref. 34], and the equilibrium solu- 
tion is compared with Blottner's and Smith's results as presented by 
Lewis [ref. 33]. 

Figure 28 shows the boundary-layer displacement thickness pre- 
dictions for the perfect and equilibrium air cases. The present method 
of solution is seen to be in excellent agreement with Adams' solution 
for the perfect gas case. The present equilibrium air solution is in 
good agreement with the results given by Keltner and Smith [ref. 35] 
but differs by a factor of two from the solution given by Blottner 
[ref. 1]. It is noted that Blottner used finite rate chemical re- 
actions and complete multi-component diffusion whereas Keltner and 
Smith used a finite-rate binary air model. 

Plots of the skin-friction coefficient and Stanton number distri- 
butions are shown in Figures 29 and 30. For the perfect gas case the 
skin friction predicted by the present method of solution is in excel- 
lent agreement with Adams' solution for S/r^ 0.4, but the solutions 
differ by up to 10% near. the stagnation point. The Stanton number 
distribution given by the present method is in good agreement with 

Adams' results for all S/r . 

n 

For the equilibrium air solution, the present method is in excel- 
lent agreement with Blottner's solution for both skin-friction and 
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Stanton number distributions 


The maximum differences between Blottner's, Keltner and Smith's, 
and the present results are less than 5%. 

7.2.3 Equilibrium and Perfect Gas Solutions for Laminar and 
Turbulent Flow Over a Spherically-Blunted Cone 
A 5° half-angle spherically-blunted cone at an altitude of 70,000 
ft. is considered in this section.* The free-stream Mach number and 
temperature were 19,564 and 226.98°K respectively. The wall tempera- 
ture was assumed to be constant at 1000°R. The problem data are pre- 
sented in Table VII, The geometry and pressure distribution are shown 
in Figures 31 and 32. The pressure distribution for this case was 
determined for perfect air, y * 1.4, using the computer program devel- 
oped by Inouye, Rakich, and Lomax [ref, 36], 

This problem is included to illustrate the solution of a weak 
adverse pressure gradient resulting from an overexpansion and recom- 
pression of the gas. Solutions are given for both laminar and fully 
developed turbulent flow for specific heat ratios, y» of 1.4 and 1.165 
and equilibrium air. 

The boundary-layer displacement thickness, skin-friction, Stanton 
number, and heat-transfer rate distributions for these solutions are 
presented in Figures 33-37. Since experimental data or other numerical 
solutions are not available for comparison, it is not possible to 
assess the accuracy of the results. However, because of the similarity 

*This case was provided by Mr. Randy Graves of the NASA Langley 
Research Center as representative of a body of unpublished free-flight 
experimental data available for this body. 
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with the hyperboloid discussed in Section 7,2,2, these results appear 
to show the correct trends, 

7,2,4 Turbulent Flow of Equilibrium Gases in Axisymmetric Nozzles 

The data for this case are given in Section 7.1.2.C of this 
Chapter, Edenfield [ref, 17] considered the nozzle discussed in this 
section for preliminary investigations leading to the design of con- 
toured nozzles for hypersonic hotshot wind tunnels with a test sec- 
tion Mach number of about 20. Edenfield used a number of theories to 
predict the downstream boundary-layer displacement thickness, but the 
available methods of prediction either failed to give a complete solu- 
tion for the nozzle or the results were found to be unacceptable. 

Attempts to solve the problem using the Elliott, Bartz, and Silver 

* 

[ref. 15] integral method failed at z/r = 1350.* This failure was 
attributed to the assumed power law total enthalpy profiles used in the 
integral method. All other attempts to predict the displacement thick- 
ness used the momentum equation only with the Crocco enthalpy distribu- 
tion or correlation formulas. 

Before discussion of the present results, a brief description of 
the physical nozzle contour and the contour used in the calculations 
is given. The physical nozzle consists of a subsonic region composed 
of a 10° conical section followed by a 5° conical section. The throat 
is cylindrical with a length of one throat diameter and the supersonic 
portion is a 5° conic section. For the calculations, the nozzle con- 
tour was smoothed as shown in Figure 17, 

' ■ » — — — ■ TO— A 

*In this section r denotes the non-dimensional throat radius. 
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It should also be noted that the total length of the nozzle was 
given as 92.087 inches (Edenfield [ref* 17]. Using a throat radius of 
0.055 inches results in a maximum value of z/r of approximately 1684. 
However, experimental pressure measurements are given in Figure 19 for 
values of z/r* up to approximately 1850. The pressure data used up- 

•ff 

stream of z/r ■ 100 corresponds to an isentropic expansion of the 
gas (nitrogen) « The transition from the isentropic expansion pressure 
data to the experimental data is particularly noticeable in the larger 
scale plots of the boundary-layer thicknesses computed using the pre- 
sent method of solution for z/ r between 100 and 400. With reference 
to the above comments, it is not clear what conclusions should be 
made in the comparison of the predicted and experimental displacement 
thickness. However, the same conditions have been used for all pre- 
diction methods and the results of the different solutions may be 
compared. 

The momentum thickness distribution corresponding to the Elliott, 
Bartz, and Silver integral method of solution and the present method 
is shown in Figure 38, The velocity or boundary-layer and displacement 
thickness are shown in Figure 39, and the predicted heat-transfer rates 
are given in Figure 40. The results presented in Figure 40 from the 
integral method were obtained by solving both the momentum and energy 
equations with a power law for the total enthalpy profiles. 

The differences in the boundary-layer thickness predictions vary 
from 10% for the velocity thickness to a factor of 2 for the momentum 
and displacement thicknesses. The heat-transfer rates predicted by the 
Integral method are from 30 to 60% lower than the rates predicted by 


71 



the present method* The results obtained with the Integral method do 
not appear to be realistic, and as noted previously the method failed 
at z/r* - 1350. The assumed power law total enthalpy profiles to- 
gether with the assumed velocity profiles in the integral method 
resulted in solutions which predicted large densities near the outer 
edge of the boundary layer in the downstream region of the nozzle. 

The density profiles predicted by the integral method, the method of 
Enkenhus and Maher [ref. 37] using the Crocco enthalpy distribution, 

■k 

and the present method are shown in Figure 41 for z/r ■ 984. As 
noted by Edenfield [ref. 17], pitot pressure measurements would detect 
the presence of these peaks in density if they exist, but the measure- 
ments made indicate that such peaks do not exist. The density profile 
predicted with the present method shows a decrease in density near the 
wall which is characteristic of boundary layers for cold wall cases. 

Figure 42 shows the boundary-layer displacement thickness pre- 
dicted by Enkenhus and Maher [ref. 37], Burke and Bird [ref. 38], Lee 
[ref. 39], Elliott, Bartz, and Silver [ref. 15], and the present 
method of solution. The present two-dimensional displacement thick- 
ness prediction is in general agreement with Burke's calculation based 
on the edge Reynolds number. The axisymmetric value of the displace- 
ment-thickness prediction is in general agreement with the solution 
of Lee. Burke's solution based on a reference Reynolds number (the 
reference Reynolds number is computed using Eckert's reference condi- 
tions) gave the best agreement with the data presented near the nozzle 
exit, but as noted at the beginning of this section, the experimental 
data presented at z/r - 1825 is outside of the nozzle dimensions. 
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The methods of Enkenhus and Maher, and Elliott, Bartz, and Silver are 
in poor agreement with the experimental data presented and the other 
solutions. 

Figure 43 shows the displacement thickness predictions of the 
Elliott, Bartz, and Silver method for both an assumed power law and the 
Crocco enthalpy distribution. For the Crocco enthalpy distribution 
only the momentum equation is solved by the integral method. Two solu- 
tions are given corresponding to different skin-friction laws using the 
method of Enkenhus and Maher, The present two-dimensional and axisym- 
metric solutions are in good agreement with the other solution proce- 
dures which use the momentum equation only. It is not clear why these 

prediction methods should be in good agreement with experiment for 
* 

z/r up to 1100 and overpredict the displacement thickness by a factor 
of two at the nozzle exit. 

To demonstrate the feasibility of obtaining solutions to problems 
with more severe conditions than the cases presented, the nozzle 
geometry, wall enthalpy distribution and reservoir density for the 
above case have been retained and the reservoir pressure and tempera- 
ture were changed to 58, 522 psia and 18,000 o R respectively. Solutions 
were obtained for perfect nitrogen, equilibrium nitrogen, and equili- 
brium air for fully developed turbulent flow. The pressure distribu- 
tion corresponds to a one-dimensional expansion of the gas and is 
shown for the three cases in Figure 44. The exit Mach numbers for per- 
fect nitrogen, equilibrium nitrogen, and equilibrium air were 21,64, 
13,69, and 10,62 respectively. 

The predicted heat-transfer rate and Stanton number distribution 
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are shown in Figures 45 and 46. The heat-transfer rates predicted for 
the equilibrium gases were from two to five times greater than the 
perfect gas solution for the region near and downstream of the nozzle 
throat. As a result of the different edge conditions , the Stanton 
number distribution shows less variation. 

7.3 Perfect Gas Solutions for Laminar and Turbulent 
Flow Over a Flat Plate With Normal Mass Injection 
Three flat-plate solutions are presented for air into air mass 
transfer and one case with no mass transfer corresponding to the exper- 
imental data referred to as run no's. 8, 11, 15, and 19 by Danberg 
[ref. 40], The free-stream Mach number for these cases was approxi- 
mately 6.3 and the temperature was approximately 59°K. The injection 

-4 

paramater, c “pv/pu, varied from 0 to 25.8 x 10 . The ratio of 

q w w e e 

the wall-to-edge static temperature was approximately 4.0. More com- 
plete data are given in Table VXIX. 

The velocity profiles obtained numerically are compared with 
Danberg’s experimental data for the four cases in Figure 47. The num- 
erical solutions correspond to fully developed turbulent flow where 
the inner eddy viscosity law is the modified Van Driest expression 
given by Cebeci [ref, 13], For y < 6, the predicted velocity is 
consistently lower than the experimental data, and except for the case 
with no mass injection, run no, 19, the velocity derivative at the 
wall given by the numerical solution is considerably less than the 
corresponding experimental value. The computed velocity profiles are 
not in good agreement with the experimental data for any of the cases 
considered having mass transfer. 

The predicted skin-friction and heat-transfer distributions are 
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shown in Figures 48 and 49. The predicted heat-transfer distribution 
is in good agreement with the experimental data. However, as would be 
expected from the poor agreement of the predicted and experimental 
velocity profiles, the agreement between the predicted and experi- 
mental skin-friction data is unsatisfactory. It is noted that the 
experimental heat-transfer data were determined indirectly by the use 
of thermocouples imbedded in the plate. The experimental skin-friction 
data were determined by two-point differentiation of the experimental 
velocity profile and are more likely to be influenced by the presence 
of the probes within the boundary layer. Tabulations of the predicted 
and experimental data for heat transfer and skin-friction are given in 
Table IX for the turbulent solution, 

7.4 Convergence Test and Computing Time Requirements 

The solutions presented above have been obtained with e^ - 0.01 
using the function test and the additional requirement that the mini- 
mum number of iterations at the new station be 2 or greater. The last 
requirement was found to be necessary to insure a smooth distribution 
of the displacement thickness when considering turbulent flow problems 
with weak pressure gradients. 

To determine the influence of upon the computing time and the 
resulting solution, a number of the above cases were computed with 
* 0,001, The flat plate Case I considered in Section 7.2.1 re- 
quired approximately ten times more computing time with *» 0.001 
than with * 0.01, The results of the two solutions differed by 
less than 2 %, Essentially the same results were obtained for other 
body shapes where the flow was turbulent and the ratio of wall to 
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stagnation enthalpy was small. The computing time was not significantly 
influenced by e^ for laminar flows with zero pressure gradient, since 
these flows are self similar. 

Representative computing time requirements for an IBM 360-65 
computer are given in the table below. The times given are for solu- 
tions where a large number of stations were solved and should be re- 
garded as maximum computing times rather than the optimum. The tabu- 
lated times for turbulent flows are for the use of the Reichardt inner 
eddy viscosity law. The solution times would be increased by approxi- 
mately 20 to 30% using the Van Driest inner law for the same number of 
stations as when using the Reichardt law. 





Execution Time 

(Min: Sec) 





Laminar 

Turbulent 

Body 

Section 

Perfect 

Equilibrium 

Perfect 

Equilibrium 

Shape 

Number 

Gas 

Gas 

Gas 

Gas 


7.2.1 
(Case 1) 

0 : 33 

- 

1:32 

2:4 

Flat Plate 

7.2.1 
(Case 2) 

0 : 32 

- 

1:9 

2:36 

Hyper- 

boloid 

7.2.2 

0 : 51 

15:1 



Blunt Cone 

7.2.3 

2 : 27 

11:22 

4:9 

19:3 


7.1.2 (B) 

- 

- 

2:26 

- 

Nozzle 

7.2.4 



7:3 

22:50 
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VIII, CONCLUSIONS AND REMARKS 


The results which have been presented demonstrate that the Crank- 
Nicolson type Implicit finite-difference scheme can be used success- 
fully to obtain solutions of the turbulent boundary-layer equations 
for nozzle flows of perfect gases or mixtures of perfect gases in chem- 
ical equilibrium. For the low Mach number perfect gas solutions con- 
sidered for nozzle flows, the present method of solution was in excel- 
lent agreement with the experimental data. For the hypervelocity 
nozzle, the present solutions were in good agreement with other numeri- 
cal solutions which used only the momentum equation and the Crocco 
enthalpy distribution. However, these solutions were not in good 
agreement with the available experimental data, but these experimental 
data did not match with the given geometry near the nozzle exit. The 
equilibrium solution for the hypervelocity nozzle was found to be 
essentially the same as that for a perfect gas at the same reservoir 
conditions. Thus, the data which were available do not appear to be 
sufficient to establish the accuracy of the numerical method. 

The results presented for equilibrium laminar flow over a hyper- 
boloid were found to be in excellent agreement with other numerical 
solutions for the predictions of heat transfer and skin friction, but 
the displacement thickness calculations differed by a factor of two or 
more among the different numerical solutions which were available for 
comparison. 

Solutions of flat-plate laminar and turbulent boundary-layer flows 
of an equilibrium gas were obtained, but these solutions were also 
essentially the same as for a perfect gas. The solutions obtained were 
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in good agreement with the experimental data. 

In all of the turbulent boundary-layer calculations considered 
without mass transfer at the wall, the solutions obtained were found 
to be essentially Independent of the expression used for the inner 
eddy viscosity law. That is the pressure gradient term appearing in 
the inner eddy viscosity law may be included or excluded without sig- 
nificantly influencing the results. 

For flat-plate turbulent flows with mass injection at the wall, 
the numerical results were in good agreement with the indirectly 
determined experimental data for heat transfer, but the predictions 
of skin friction and velocity profiles, determined from direct measure- 
ments, were not in satisfactory agreement with the experimental data. 
Thus, these results are also inconclusive. 

The present method of solution has been demonstrated to be accu- 
rate for both internal and external flows of perfect gases for both 
laminar and turbulent flows for the case of no mass injection. More 
experimental data are needed to verify the method for the case of mass 
transfer. 

For the equilibrium chemistry calculations, the present method of 
solution for the laminar flow of the gas is in good agreement with 
other numerical solutions. For turbulent flow of an equilibrium gas 
the present method of solution was found to be stable and the solutions 
show the expected trends, but the experimental data which were avail- 
able were in a range of pressure and/or temperature which showed little 
departure from the perfect gas condition. With the exception of the 
Elliott, Bartz, and Silver [ref# 15] integral method of solution for 
equilibrium chemistry turbulent boundary-layer flows in nozzles, other 
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numerical solutions were not available for comparison. For the case 
of a hypervelocity nozzle, the integral method failed to converge for 
a stagnation pressure and temperature of 11,511 psia and 3000°K. 
Whereas, the present method was found to be stable for nozzle calcula- 
tions with reservoir pressures and temperatures up to 58,522 psia and 
10,000°K respectively. 
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IX. APPENDIX A 


Thermodynamic and Transport Properties and one Dimensional 
Expansion of a Reacting Gas Mixture In Chemical Equilibrium 

In order to solve the conservation equations for laminar or tur- 
bulent boundary-layer flow of an arbitrary mixture of perfect gases in 
thermodynamic equilibrium, it is necessary to specify the boundary- 
layer edge conditions and to provide a method for determining the 
local thermodynamic and transport properties for the gas mixture con- 
sidered. 

The thermodynamic and transport properties and the inviscid ex- 
pansion data which have been used in the present solution of the 
boundary-layer equations have been obtained by modifying the basic 
computer program developed by Lordi, Mates, and Moselle [ref. 20], 
which was developed to determine the inviscid expansion of arbitrary 
gas mixtures assuming either frozen chemistry, equilibrium chemistry, 
or finite-rate chemistry. Only the equilibrium chemistry option is 
considered in the sections below. 

A.l Equilibrium Gas Properties 

A brief description of the procedure for determining the equili- 
brium chemistry solution is included in this appendix. The notation 
is the same as that employed by Lordi, Mates and Moselle, and is in- 
cluded in a separate list of symbols presented at the end of this appen- 
dix. The analysis presented below is essentially the same as that 
given in Sections 2.1, 2.2, and 4,1 of their report. 
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A. 1.1 Thermodynamic and Transport Properties 


The enthalpy of a given gas mixture in thermodynamic equilibrium 
is a function of the internal energies of the species and the chemical 
composition of the mixture, which are related to the local thermo- 
dynamic properties, pressure, density, and temperature. The composi- 
tion of a specified gas mixture at a given equilibrium state is 
determined for c chemical elements and s chemical species as follows. 
The chemical equation for the ith species may be expressed as 

— M^(ct^, a ±2** * * * a ic^ * ^ *" 2, ....,8 (Al) 

The species to be included in the mixture are specified through the 
matrix ct^. 

The equilibrium formation reactions for the species in terms of 
the elements may be expressed by .(s - c) linearly independent rela- 
tions of the form 

c 

M. — £ ct . . M , , i — 1, 2,....s (A2) 

1 j-l 30 J 

where the rank of the a . matrix is c, Lordi, Mates and Moselle note 

XJ 

that in some cases it is necessary to choose species other than the 

chemical elements as components. The formation reactions for the 

(s - c) dependent species are expressed in terms of the c components as 

c _ 

M. * £ v. . M., i * c + 1, c + 2,,,..,s (A3) 

1 j-l 2 2 

The number of gram-atoms of each chemical element, Q^» contained 
in the mixture must be given in addition to the chemical species and 


81 



chemical elements# The total number of gram-moles of the components 
which are present in the mixture is given by 

X “jk ij ■ \ (M) 

This relation specifies the element composition in terms of the com- 
ponents rather than the chemical elements. 

Enforcing mass conservation in each of the formation reactions and 
a global mass conservation gives the mole fractions of the components 
as 

s _ 

X. “ q. - l tv.. - q. (v - 1)] X , j - c + 1, c + 2, . ,s (A5) 
i i i-c+1 ^ i 1 1 


where 






The dependent species concentrations are related to the component con- 
centrations through the equilibrium constants for the formation re- 
actions, equation (A3), by 


Y 


P' 


(v ± -l) 


c 

n 


ij 


j-1 


c+1, c + 2, < 


(A6) 


where 



exp 


o o 
U# P 


L v ij T T 


j-1 


(A7) 
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The concentrations of the independent species or components are 
determined by solving equation (A5) and using a Newton-Raphson itera- 
tion procedure to solve equation (A4) . The molecular weight at the 
specified pressure and temperature is calculated from 

s 

M - l X, M (A8) 

j-1 3 3 

The density is then determined from the equation of state for a mix- 
ture of perfect gases 

R 

P' “ P' T* (A9) 


and the specific enthalpy is given by 


H - l X h 
j-1 3 3 


(A10) 


For a mixture of perfect gases, the molar enthalpy, h^ , chemical poten- 
tial, at one atmosphere of pressure, and molar entropy, s°, are 
functions of temperature only. These data are specified by either a 
simple-harmonic-oscillator model, or by polynomial curve fits to more 
accurate calculations of the species properties. The latter method is 
called the thermo-fit method. In both methods, the species are assumed 
to be vibrationally and electronically excited, in equilibrium with the 
local translational temperature. The translational and rotational 
degrees of freedom are assumed to be fully excited. 

For the harmonic-oscillator model, the specific enthalpy for a 
monatomic or diatomic species is expressed as 
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5 + 2(n_j - 1) 
2 


+ (n 


j 


1) 



I e g exp[- e /T] 

■ 1 - , *■!- j - -1 - ■ 

exp(6 v /T)-l Lj 

j T l g [- e /T] 

£-1 j j 


and the chemical potential at standard pressure is given by 



{ 5 + 2(n - 1) 

j a j + 2 ln T 


expt9 v /T] 

+ - l) ln expl'cV'/T) - 1] 

3 


+ In 


£-1 



exp[- z z /T] V 


where 


5 + 2(n - 1) 

“j * bj + 2 ln T 'o 


and 


b 


j 



2ir M' k 

■ Mi la J l ■ II 


+ In 



(n . - 1) ln 0' 

j r j 


(All) 


(A12) 


(A13) 


(A14) 
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For the thermo-fit method, these data are given by 


o , o 

“J ‘ h 3. 


- a.(l - In T') - 


c d 6 

bj T' (T*) 2 - -1 (T*) 3 - -1 (T') 4 - k 


and by differentiation of equation (A15) 


(A15) 


h j “ h j 


« a + bj T* + Cj (T*) 2 + dj (T 1 ) 3 + * 3 (T') 4 (A16) 


For both methods is determined from 


o Vi 

3 » mil l m m 0l 


(A17) 


and the specific heat of the individual species is determined by dif- 
ferentiation of equation (A12) or equation (A16). 

The local thermodynamic and transport properties which must be 


specified for use in the boundary-layer computer program ares A T', 

c„ 


'V/ 

H 


log 


10 


H' 


EGAS 


P * Log 10 (P'), p - Log 1Q (p /p s ), p , Cp 


RGAS 


, Pr. 


Tables of these data are generated for the desired range of 
temperatures in decreasing T 1 at constant pressure. The procedure is 

<\j 

repeated for decreasing in P to span the desired range in pressure. 
Conversion of the data to the required form is as follows 


H ■ log 


fHT* M 1 
s 


10 


M 


log 10 (H»/RGAS) 
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'V 

p 


log 


p' (2.205 x 1(T 3 ) 
o 


10 (3.531 x 10” 5 ) (32.176) (p ’) 

s 


f3H 

" M 


& 

10 /log 10 (e) 


The derivative is computed using a five point Lagrangian differentia- 
tion polynomial. 

The Prandtl number is obtained by interpolation in the tables of 
data given by Hansen [ref. 41] for air and Ahtye and Peng [ref. 42] for 
nitrogen. 

The viscosity of the gas mixture is computed using Wilke's [ref, 
43] semi-empirical formula 


* 

U - 


ISS 


l 

i-1 


Y i u i 


ISS 


U-1 * 1J l 


where 



- ' 

-1/2 


*1 

1/2 

’ ' 

1/4 

1 + -=• 


1 + 

fc 


!i 





v i. 





(A18) 


(A19) 


Species viscosities f are approximated by curve fit data of the 

form 
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* C (A In T + B ) 

■ e T gm/cm-Sec (A20) 

•k 

The constants A^> and are obtained from curve fitting the 
data given by use of the method of Yun and Mason [ref. 44]. The data 
which have been used in the boundary-layer calculations were given by 
Blottner [ref. 7]. 

The frozen values of thermal conductivity and Prandtl number, and 
the compressibility factor are also computed, but these quantities were 
not used. 

A description of the modifications made to the computer program 
of Lordi, Mates and Moselle (for computing the thermodynamic properties) 
and the mixture transport property calculations are given by Miner, 
Anderson and Lewis [ref. 28] with listings of the added or modified 
subroutines, and a description of the necessary input data is given. 

A. 1.2 Boundary-Layer Edge Conditions for a 
Gas Mixture in Thermodynamic Equilibrium 

The boundary-layer edge conditions for the equilibrium gas case 
are obtained from the isentropic expansion of the gas from a given 
stagnation or reservoir state. To obtain the inviscid expansion of 
the gas it is necessary to conserve momentum and energy in addition to 
the requirements given in Section A, 1.1* The procedure for determining 
the inviscid expansion discussed below' is essentially that given in 
Section 2.2 of Lordi, Mates and Moselle [ref. 20]. 

For a quasi-one-dlmensional inviscid flow of a gas the conserva- 
tion of mass, momentum, and energy is expressed as 
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puA ■ M ■ constant 


(A21) 


udu + — dp «* 0 (A22) 

P 

2 

H + “ H q “ constant (A23) 

For a chemically reacting gas mixture the entropy is expressed as 




s 

l 




- I X In X - In p' 
j-1 3 3 


constant 


(A24) 


for an isentropic path. The governing equation which must be satisfied 
by the c components and the pressure is 


Fj *= 0, j - 1, 2,...., c + 1 

where from (A5) 

s __ _ 

F. " q. - l [v . - q. (v. - 1)] X. - X., j « 1, 2...C (A25) 
3 3 i»c+l 13 3 1 1 3 


and from (A24) 


F 


c+1 


8 

I X.(s° - M. S » - In X.) - In p* 
i-1 1 1 1 ° ± 


(A26) 


The expansion is obtained by taking successive temperature steps 
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from the reservoir value and solving equations (A25) and (A26) by a 
Newton-Raphson method to determine the composition and pressure at the 
specified temperature and entropy. 

To express the solution in terms of area ratio, the critical mass 
flow is first determined by computing the maximum value of pu. This 
maximum value with the corresponding temperature determines the throat 
conditions for nozzle flows. The solution is then restarted from the 
reservoir and the area ratio is determined from equation (A21) , At 
each point in the expansion, the density and Mach number are computed 
from 


R 



M 


and 


M - u/[(P - P b )/(p - P b )] 1/2 


where the subscript b refers to the values at the previous step. 
The area ratio is computed using 


(pu) 


max 


pu 

The velocity, pressure, and enthalpy are converted to the form 
required by the boundary layer program as follows 

u* - u(R (778.158) (1.8) (32.176) T*/M) 1/2 
o o 


* ■ ** 1 o ( **P 
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a - log 10 (H r M s /M) 


Lordi, Mates and Moselle's computer program [ref. 20] was modified and 
is described by Miner, Anderson, and Lewis [ref. 28] where the addi- 
tional input data are given. 
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LIST OF SYMBOLS 


constant in thermo-fit and harmonic oscillator expressions 
for chemical potential, 

ratio of local cross-sectional area to minimum cross- 
sectional area, 

coefficient of T* in thermo-fit expression for species 
enthalpy, or defined by Eq.(Al4) for harmonic-oscillator 
description, 

number of elements in mixture. 

2 

coefficient of (T*) in thermo-fit expression for species 
enthalpy, or mass fraction. 

2 2 

specific heat of mixture (ft /sec -°K), 

* 

dimensionless specific heat (c p /RGAS). 

3 

coefficient of (T') in thermo-fit expression for species 
enthalpy , 

coefficient of (T')^ in thermo-fit expression for enthalpy, 
degeneracy of the 1th electronic state of the jth species. 


Flank's constant, 

molar enthalpy of the jth species (h_j « h^/R Q T^) 


specific enthalpy of mixture 


H » 


H* M ' 
o 

R T* 
o o 


log 


fHT* M 1 
s 


10 


M 


or log 10 (H'/RGAS) 


Boltzmann's constant 




M. 


constant in thermo-fit expression for the chemical 
potential of the jth species. 

eqjilibrium constant for the 1th reaction based on partial 
pressure 

mass of the jth particle (gm). 


M 

ft 

M 


*3 

P 

P 




RGAS 

s 

S j 

T 

* 

u 


molecular weight (gms/mole) . 

Mach number 

molecular weight of the gas mixture at standard atmos- 
pheric conditions, 
number of atoms in the jth species. 

pressure (P * P'/P^), P’ in atmospheres. 

io8io<pp;> 

mole fraction of jth component when only independent 

species are present in mixture. 

number of gram atoms of kth element in mixture. 

universal gas constant 1.98647 

gm-mole K 

2 2 

gas constant (ft /sec -°K) 

number of species in mixture 

molar entropy of jth species (s. “ s'/R ) 

3 jo 

temperature (T’/T^); T’ in °K 
velocity (ft/sec) 


velocity ju ** u ' / J T VM q 
mole-fraction of jth species 


92 



matrix specifying the elements and species in the chemical 
model of the mixture 

mass concentration of the ith species, lb-mole/lb 


energy of the £th electronic state of the jth species 


“j 


'I/W 


characteristic relational temperature of the jth species 


characteristic vibrational temperature of the jth species 


(0 « 0* /T 1 ) 

V j "J ° 


1 


V j " R T* 


o o 


chemical potential of the jth species 
species viscosity (gm/cm-sec) 

stoichiometric coefficients of the equilibrium formation 
reactions 


j 


ij 


density (p - p'/p^) 
3 

density (slugs/ft ) 
logi 0 (P*/p*) 


Subscripts 
reservoir condition 
standard atmospheric condition 
pertaining to the ith reaction 
pertaining to the jth species 
pertaining to the kth element 
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pertaining to the 1th electronic level 


Superscripts 
dimensional quantity 
refers to standard pressure condition 
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TABLE I 


Test Conditions for Cole's Flat Plate Experiments 


Case 

P G (cm Hg) 

T q (°R) 

u (ft/sec) M 

CO CO 

Re xlO ^ /ft 
00 

T 

w 

20 

103 

561 

2219 

3.701 

0,1642 

Adiabatic 

26 

131 

568 

1972 

2.578 

0*387 

Adiabatic 

62 

302 

563 

2327 

4.544 

0.322 

Adiabatic 


-99 



TABLE II 


Elliott* Bartz. and Silver 
P q “ 300 psia 

T » 4500 °R 
o 

Re “ 0.2 x 10^/in, 

CO 

M «* 2.685 

co 

u • 6231 ft/sec, 

CO 

r - 1.2 


Teat Case Nozzle Data 
P * 11.54 psia 

T - 2614 °R 

09 

« 3,7 x 10 ^ slugs/ft^ 

T - 1125 °R 
w 

Pr - 0.83 
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TABLE III 



M i 

CD 

Y - 


Test Conditions for NASA-Lewia Nozzle - Case 2a 
1 299 « 3 psia * 0.56 psla 

■ 970 °R T - 161 °R 

CO 

“ 0*588 x 10^/in, « 2.9 x 10 ^ slugs/ft^ 

1 5.01 u * 3118 ft/sec 

00 


1.4 


T variable 
w 
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TABLE IV 


Test Conditions for AEDC Hotshot-Nozzle 
Perfect Gas C^) 

P q ** 11511.57 psia P^ ■ 4.34 x .10 4 psia 

T - 5400 °R T - 78.87 °R 

o ® 

Re “ 5.3 x 10 4 /in. p ■ 4.46 x 10 ^ slugs/ft^ 

CO CO 

M •» 18.37 u •• 8132 ft/sec 

CO CO 

y ■ 1.4 T , variable 

w 


Equilibrium 

P - 11511.57 psia P 

o 00 

T - 5400 °R T 

O 00 

Re - 3.0 x 10 4 /in, P 

CO 00 


-3 

4.33 x 10 psia 
108 °R 

3.41 x 10 6 slugs/ft^ 


M * 17.24 u « 8732 ft/sec 

CO CO 

T . variable 

w* 
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TABLE V 


Test Conditions for Hironimus’ Flat Plate Experiments 


Case 

P « 8498 psia 
o r 

T - 4913 °R 
o 

Re ■ 4.46 x 10^/in 

00 

M - 7.391 
00 

T - 533 °R 
w 

P « 14670 psia 
o 

T - 3691 °R 
o 

Re - 1.33 x 10 6 /in 

CO 

M ■ 7.58 

CD 

T - 355 °R 
v 


P ro « 1.53 psia 
T - 546 °R 

CO 

p * 2,49 x lCf 4 slugs/ft^ 

CO 

u « 8157 ft/sec 
00 

P ra - 2.46 psia 

T - 339 °R 
00 

p ro ■ 6.1 x 10 -4 slugs/ft^ 
u • 6843 ft/sec 

OO 
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TABLE VI 


Flight Conditions 

P* " 6.035 atms 
o 

T' - 6996 °K 
o 

Re - 2157943/ft 

CO 

M - 20.178 

CO 

u ■ 20,000 ft/sec 

CO 


for AGARD Test Case A 

P - 1.0997 x 10 
00 

T - 226.98 °R 

OO 

ALT - 100,000 ft 
r «* 1 in 

IX 

T - 1400 °K 
w 


atms 
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TABLE VII 


Flight Conditions for Sphere-Cone Sample Case 


P'/P - 513.994 

O 03 

T'/T - 33,065 
o « 

p'/p. - 11.783 

O 00 

M - 19.564 

CO 

Re - 8.898 x 10 6 /ft 
00 

u “ 19,000 ft/sec 

CO 9 


P - 4.429 X 10 

CO 

T - 217.9 °K 

CO 

p - 5.555 x 10‘ 

oo 

Alt - 70,000 ft 

T - 1000 °R 
w 

y’ - 1.1655 
o 


atms 


AMAGAT 


105 


TABLE VIII 


Test Conditions for Flat Plate Flows with Mass Injection 


Run 

No. 

x, mm 

Re x x 10”^ 

M 

e 

T /T 
w e 

T °K 
e 

P v 

w _W . ‘ 

p u 
e e 

5 

377.8 

2,925 

6.38 

4.102 

60.4 

9.23 

6 

428.6 

3.446 

6.49 

4.132 

58.1 

9.07 

7 

479.4 

4.160 

6.60 

4.263 

56.3 

8.85 

8 

530.2 

4.261 

6.48 

4.167 

58.3 

9.21 

9 

377.8 

3.249 

6.39 

4.102 

59.5 

16.86 

10 

479.4 

3.950 

6. .42 

3.991 

59.1 

16.24 

11 

530.2 

4.264 

6.37 

3.808 

60.0 

17.81 

12 

377.8 

3.185 

6.28 

3.864 

61.9 

24.91 

13 

428.6 

3.619 

6.20 

3.640 

62.4 

24.52 

14 

479.4 

3.945 

6.22 

3.759 

62.5 

24.90 

15 

530.2 

4,298 

6.22 

3.662 

62.5 

25.80 

16 

377.8 

3.057 

6.52 

4.418 

57.6 

0.0 

17 

428.6 

3,163 

6.32 

4.171 

61.0 

0.0 

18 

479.4 

3.679 

6.44 

4.346 

59.0 

0.0 

19 

530.2 

4.260 

6.45 

4.297 

58.6 

0.0 
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TABLE IX 


Comparison of Numerical Predictions and Experimental Data for Flat Plate Flows with Mass Injection 


Danberg 
Run No. 

Exp. or 
Calc. 

x(mm) 

Re /mm 

p v /p u 
w w e e 

C f xl0* 

C.xlO 

h 

0 (mm) 

6 (mm) 

6 (mm) 

Re 0 

6*/0 


exp. 




12.54 

6.33 

0.376 

5.68 

10.0 

3023 

15.11 

19 


530.2 

8034.7 

0.0 









calc. 




11.95 

7.77 

0.409 

5.75 

10,1 

3284 

14.06 


exp. 




8.39 

3.17 

0.658 

7.86 

14,1 

5288 

11.94 

8 

calc. 

530.2 
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Figure li Boundary Layer Coordinate System 
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Figure 4: Velocity Profiles for a Flat Plate-Law of the Wall 
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Figure 7: Nozzle Geometry-Elliott, Bartz, and Silver Test Case 
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Figure 8t Eddy Viscosity Profiles at Nozzle Throat-Ell iott, Bartz, 
and Silver Test Case 
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Figure 9: Heat Transfer Distribution Along Nozzle Wall-Elliott, Bartz, and Silver Test Case 
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Figure 10» Boundary Layer Thickness Distribution for a Nozzle-Elliott, Bartz, and Silver 
Test Case 
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Figure 12 t Momentum Thickness Distribution for a Nozzle-Ell iott, Bartz, and Silver Test Case 
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Figure 13s Geometry for the NASA-Lewis Nozzle 


6 

z, in. 











0 2 4 6 8 10 12 

z . in. 


Figure 15 t Wall Temperature Distribution for the NASA-Lewis Nozzle 
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Figure 16 » Heat Transfer Distribution for the NASA-Lewis Nozzle 





Figure 17s Geometry for the AEDC Hotshot-Nozzle 




Figure 18 t Wall Enthalpy Distribution for the AEDC Hotshot-Nozzle 



Figure 19 l Pressure Distribution for the AEDC Hotshot-Nozzle 
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Figure 21 t Heat Transfer Distribution for the AEDC Hotshot-Nozzle 
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Figure 23 t Number of Iterations Required to Obtain a Converged Solution 



Figure 24* Skin Friction Distribution for a Flat Plata-Hironimus * Data Case I 






for a Flat Plate-Hironimus * Data Case II 
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Figure 26: Geometry for a 10° Half-Angle Hyperboloid-AGARD Case A 



Figure 27 t Pressure Distribution for a 10° Half-Angle Hyperboloid-AGARD Case A 
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Figure 28 i Displacement Thickness Distribution for a 10° Half-Angle Hyperboloid-AGARD Case A 





Figure 29 t Skin Friction Distribution for a 10° Half-Angle Hyperboloid-AGARD Case A 




Figure 30s Stanton Number Distribution for a 10° Half-Angle 
Hyperboloid-AGARD Case A 
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Figure 31 l Geometry for a 5° Half-Angle Spherically-Blunted Cone 
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Figure 34 J Displacement Thickness Distribution for Turbulent Flow 
Over a 5° Half-Angle Spherically-Blunted Cone 
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Figure 35 t Skin Friction Distribution for a 5° Half-Angle Spherically-Blunted Cone 
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Figure 37: Heat Transfer Distribution for a 5° Half-Angle Spherically-Blunted Cone 



Figure 38: Momentum Thickness Distribution 




z/r K 


Figure 39 i Boundary Layer and Displacement Thickness Distribution for the AEDC Hotshot-Nozzle 
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Figure 40 j Heat Transfer Distribution for the AEDC Hotshot-Nozzle 
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Figure 41? Density Profiles for the AEDC Hotshot-Nozzle 




Figure 42$ Comparison of the Predicted Displacement Thickness Distribution with Other 
Methods Using Both Momentum and Energy Equations-AEDC Hotshot-Nozzle 
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Figure 43 1 Comparison of the Predicted Displacement Thickness Distribution with Other Methods 
Which Use the Momentum Equation Only-AEDC Hotshot-Nozzle 
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Figure 44 1 


Pressure Distribution for the AEDC Nozzle Configuration-P 
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Figure 45: Heat Transfer Distribution for the AEDC Hotshot-Nozzle 

Configuration-P - 58,522 PSIA; T - 10,000°K 
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Figure 47s Velocity Profiles for a Flat Plate With Mass Injection- 
Danberg's Data 


154 


D ANBERG EXPERIMENTAL DATA, NOMINAL VALUES 
RUN NUMBERS I 16-19 5-8 9-1 I. 12-1 

SYMBOL ■ • A ♦ 


16-19 

5-8 

9-11 

. 12-15 

■ 

• 

▲ 

♦ 

6.4 

6.5 

6.4 . 

6.24 

4.31 

4.17 

3.94 

3.73 

0.0 

9.09 

17.0 

25.0 

CALCULATIONS 


yt ft 4.31 4.17 3.94 

u e x)Q 4 j 0.0 9.09 1 17.0 

q CALCULATIO 

h= /W — laminar 


c h I0 J 


FULLY TURBULENT 


RUN 19 




RUN 8 


l\ /RUN 15 
D(/RUN II 
UV/RUN 8 


7^ run ii 
>^RUN 19 

^ RUN 15 


Figure 48: Stanton Number Distribution for a Flat Plate -With Mass Injection-Danberg’s Data 
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Figure 49s Skin Friction Distribution for a Flat Plate With Mass Injection-Danberg’s Data 



